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One of the most promising proposals for engineering topological superconductivity and Majo- 
rana fermions employs a spin-orbit coupled nanowire subjected to a magnetic field and proximate 
to an s-wave superconductor. When only part of the wire's length contacts to the superconduc- 
tor, the remaining conducting portion serves as a natural lead that can be used to probe these 
Majorana modes via tunneling. The enhanced role of interactions in one dimension dictates that 
this configuration should be viewed as a superconductor-Luttinger liquid junction. We investigate 
such junctions between both helical and spinful Luttinger liquids, and topological as well as non- 
topological superconductors. We determine the phase diagram for each case and show that universal 
low-energy transport in these systems is governed by fixed points describing either perfect normal 
reflection or perfect Andreev reflection. In addition to capturing (in some instances) the famil- 
iar Majorana-mediated 'zero-bias anomaly' in a new framework, we show that interactions yield 
dramatic consequences in certain regimes. Indeed, we establish that strong repulsion removes this 
conductance anomaly altogether while strong attraction produces dynamically generated effective 
Majorana modes even in a junction with a trivial superconductor. Interactions further lead to 
striking signatures in the local density of states and the line-shape of the conductance peak at fi- 
nite voltage, and also are essential for establishing smoking-gun transport signatures of Majorana 
fermions in spinful Luttinger liquid junctions. 

PACS numbers: 73.21.Hb, 71.10.Pm, 74.78. Fk 
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I. INTRODUCTION 

Topological phases display many 'exact' features that 
are remarkably robust to variations in microscopic re- 
alization and imperfections. One of the most exotic is 
the possibility of emergent excitations known as non- 
Abelian anyons that host zero-energy 'internal' degrees 
of freedom in an intrinsically non-local way. Informa- 
tion encoded in such anyons is then naturally protected 
from decoherence and can furthermore be manipulated 
using braiding operations, forming the basis of attractive 
quantum computing platformsP Recently attention has 
focused on the so-called Majorana fermion, which binds 
to a specific kind of non-Abelian anyon originally pre- 
dicted to occur in the v — 5/2 fractional quantum Hall 
state^J. Among the many current proposals for generat- 
ing Majorana fermioniP"21 ( se e [718] for a review), one 
attractive class involves spinless one-dimensional (ID) 
topological superconductors whose hallmark is the ex- 
istence of Majorana zero-modes localized at the ends of 
the system.^ This phase can be engineered in a variety 
of settings, including 2D topological insulator edges 10 , 
spin-orbit-coupled nanowire d 11 * 12 !, 3D topological insula- 
tor nanoribbon&LS, and cold atomic gases 14 '. Numerous 
quantum computation protocols based on manipula ting 
Majoranas in ID systems have now been proposeoP^lES] 

The first step towards realizing such applications is of 
course the conclusive experimental identification of Majo- 



rana fermions. One appealing detection method involves 
transport. In particular, several studies predict that tun- 
neling electrons onto a Majorana mode gives rise to a 
zero-bias conductance anomalypsHM] Here we revisit this 
problem from a new perspective based on renormaliza- 
tion group methods similar to those of Rcfs. [36-38 , and 
especially [39 61 . A major virtue of this approach is that 
it allows one to extract universal tunneling signatures of 
these modes even when strong interactions are present. 
Indeed, we will (in some cases) recover in a very gen- 
eral way previous results based on specific model calcu- 
lations, and also identify new regimes where interactions 
lead to dramatic and very surprising consequences. Fur- 
thermore, our approach provides an elegant means of ad- 
dressing the fate of localized Majorana zero-modes when 
coupled to gapless degrees of freedom. 

The theoretical technology developed here is widely 
applicable to superconducting Majorana platforms. We 
will, however, mainly focus on spin-orbit-coupled ID 
systems such as a semiconducting nanowire subjected 
to a magnetic field, in the experimentally accessible 
geometry where half of the wire couples to an s-wave 
superconductor while the other half remains gapless; see, 
for example, Fig. [IJa) and Ref. I35|401 By gating one can 
independently tune the left and right halves between 
a 'helical' regime — with only one active channel at low 
energies — and a 'spinful' regime where two chann els play 
a role (multi-channel regimes are also accessibl e 1 34 * 41 * ^ 3 -* 
but will not be considered here). Crucially, in cither 





Helical wire 


Spinful wire (<j CT = 1) 


Topological SC 


W,9>h 


r _ J °. »p < 1 


Non Topological SC 


r J 0, g<2 
\^9>2. 


G - { °' 5p < ! 
0<G<^,<; P = 1 



TABLE I: The four junction archetypes studied in this paper. The columns denote the two possibilities for the one dimensional 
leads, and the rows differentiate between superconductors with and without Majorana zero mode bound states. The Luttinger 
parameters g and g p describe the interactions in the lead; g < 1 and g > 1 correspond to repulsive and attractive interactions 
respectively. Note that for the regime of weak repulsive interactions expected in many solid state implementations, a zero bias 
conductance G = 2e 2 jh is a robust signature of Majorana zero mode bound states in the superconductor. 



limit the gapless half generically forms a Luttinger 
liquid. The nature of the superconducting state in 
the other half of the wire depends strongly on the 
number of channels: in the helical case a topological 
phase supporting Majorana zero-modes emerges, while 
in the spinful regime a trivial gapped state appears 
instead IH112I This system therefore admits four natural 
superconductor/Luttinger liquid junction archetypes: 
the gapless part of the wire (which we view as a lead) 
can be either helical or spinful and likewise the supercon- 
ducting region can be topological or trivial. We stress 
that it is essential to understand the transport properties 
of both the topological and trivial junctions to establish 
unambiguous transport signatures of Majorana modes. 
For instance, if the conductance can behave similarly 
with or without the presence of a Majorana then clearly 
this would be less than a 'smoking- gun' detection scheme. 



To introduce the basic philosophy underlying our ap- 
proach it is useful to first imagine physically cutting the 
wire such that the Luttinger liquid and superconducting 
regions decouple entirely. Transfer of electrons across the 
junction formed by the two subsystems is then prohibited 
for trivial reasons. In renormalization group language, 
the system's low-energy behavior here is described by a 
fixed point theory characterized by a vanishing conduc- 
tance. We refer to this as a 'perfect normal reflection' 
fixed point since in this case electrons incident on the su- 
perconductor undergo normal reflection with unit prob- 
ability at the junction. Suppose now that the Luttinger 
liquid and superconductor 'reconnect', and one incorpo- 
rates arbitrary symmetry-allowed couplings between the 
two. Our objective is to then address questions such as 
the following: What is the fate of the perfect normal re- 
flection fixed point in this case? If it is unstable, what 
couplings provide the leading instability and to which 
fixed point do they ultimately drive the system? What 
are the properties of such putative fixed points? And 



what are the implications for transport experiments? 

Let us now highlight our main results (partially sum- 
marized in Table I) for the four cases that we analyze, be- 
ginning with the topological superconductor/helical Lut- 
tinger liquid junction. In this case the gapless region is 
characterized by a Luttinger parameter g where 5 = 1 
represents the free-fermion limit while g < 1 and g > 1 
respectively correspond to repulsive and attractive inter- 
actions. When the wire is cut as described above the 
topological superconductor supports a single localized 
Majorana zero-mode at the junctionPSProvided the Lut- 
tinger parameter falls in the range g > 1/2 tunneling elec- 
trons onto this mode constitutes a relevant perturbation 
that destabilizes the perfect normal reflection fixed point. 
We demonstrate that the Majorana zero-mode then de- 
localizes completely into the Luttinger liquid and drives 
the system to a fixed point describing perfect Andreev 
reflection at the junction. This perfect Andreev reflec- 
tion fixed point is characterized by the familiar quantized 
zero-bias conductance G = 2e 2 /h at zero temperature 
T = that has been captured by numerous studies in 
the free-fermion limitP^"- 30 * 32 * 33 * 35 ^ In addition to the uni- 
versal value of 2e 2 /h for the conductance at zero bias, the 
topological superconductor/helical Luttinger liquid junc- 
tion also exhibits a universal form for the finite bias con- 
ductance curve. The form of this curve can be computed 
in perturbation theory in certain voltage regimes, allow- 
ing one to extract the value of the Luttinger parameter 
g. For g < 1/2, however, (which is potentially applica- 
ble to carbon nanotube-based Majorana proposals^H^) 
coupling to the zero-mode at the junction is irrelevant at 
the perfect normal reflection fixed point, which is then 
stable. Consequently the zero-bias conductance vanishes 
at T = 0. This does not, however, imply that the Ma- 
jorana mode remains localized — we show using scaling 
that in this limit the probability density associated with 
the Majorana mode decays into the Luttinger liquid as a 
power law. 



When a helical Luttinger liquid impinges instead on 
a trivial superconductor, the perfect normal reflection 
fixed point is stable, resulting in a vanishing zero-bias 
conductance, for any g < 2. Taken together, the above 
results imply that electron transport from a helical wire 
(with not-too-strong repulsive interactions) onto a su- 
perconductor is unmistakably different in the topological 
and trivial cases and thus indeed provides an unambigu- 
ous way of identifying Majorana zero-modes. The g > 2 
regime may be realizable in cold atoms experiments^ and 
exhibits fascinating properties even from a purely theo- 
retical perspective. Here Cooper-pair tunneling from the 
Luttinger liquid to the superconductor is relevant, gener- 
ating a flow to the perfect Andreev reflection fixed point 
with quantized conductance. We show that this results 
from the dynamical generation of a pair of asymptoti- 
cally decoupled Majorana modes which mediate the per- 
fect Andreev reflection, ft is remarkable that any trace 
of Majorana physics appears in a junction with an ordi- 
nary superconductor; this would certainly be interesting 
to explore further in numerical simulations. 

Perhaps more than any other case, superconduc- 
tor /spinful Luttinger liquid junctions highlight the limi- 
tations of model-specific microscopic calculations and the 
utility of our approach based on universal low-energy 
physics. When a spinful Luttinger liquid couples to a 
topological superconductor, the leading perturbation at 
the perfect normal reflection fixed point again involves 
electron tunneling onto the Majorana zero-mode at the 
junction. Its relevance follows from the charge- and spin- 
sector Luttinger parameters g p and g„. With unbroken 
SU(2) spin symmetry in the Luttinger liquid one has 
g a — I, and in this case coupling to the Majorana is 
relevant for all g p > 1/3. Because only one of the two 
spin channels can hybridize with the Majorana mode, the 
junction flows to a novel fixed point corresponding to per- 
fect Andreev reflection for one species and perfect normal 
reflection for the other. This fixed point is in fact robust 
with respect to some perturbations which break 517(2) 
in the bulk, i.e. correspond to g a ^ I, such as spin or- 
bit coupling and magnetic field. Thus, in the topological 
case, a spinful Luttinger liquid junction again exhibits ro- 
bust G = 2e 2 /h conductance quantization down to fairly 
strong repulsive interactions with g p = I /3. 

An analysis of the non-topological superconduc- 
tor/spinful Luttinger liquid case in the free-fermion 
limit, however, yields a non-universal zero-bias conduc- 
tance ranging anywhere from to 4e 2 /h depending on 
parameters — potentially making it difficult to differen- 
tiate from the signal originating from the Majorana in 
the topological junction. This apparent non-universality 
originates from the fact that the leading perturbations 
to the perfect normal reflection fixed point in the free 
fermion case are exactly marginal. Fortunately we find 
that arbitrarily weak repulsive interactions are sufficient 
to stabilize the perfect normal reflection fixed point. 
Since weak repulsive interactions are generic in many 
physical realizations of such Luttinger liquids, the non- 



universal zero-bias conductance calculated for the free 
fermion limit will eventually renormalize to zero at suf- 
ficiently long length scales and low energies. We thus 
conclude that a quantited zero-bias G = 2e 2 /h conduc- 
tance indeed serves as a definitive fingerprint of a Majo- 
rana mode at the junction, for both spinless and spinful 
Luttinger liquids. 

The remainder of the paper analyzes each of the four 
cases in detail. We begin in Sec. |ll] with the heli- 
cal/topological junction, which we first solve in the free 
fermion limit g = 1 and then in the general case us- 
ing bosonization. In particular, we extract the conduc- 
tance from a boundary action obtained by integrating 
out the bulk of the Luttinger liquid, and use duality 
to gain insight into the nature of the two fixed points. 
In Sec. |ni we apply similar methods to analyze the 



helical/non-topological junction. Sections IV andfvltreat 
the spinful/non-topological and spinful/topological junc- 
tions. In the discussion (Sec. VI I we examine the phys- 
ical consequences of our work in more detail. In partic- 
ular, we show that the line-shape of the finite bias con- 
ductance curve contains information about the Luttinger 
parameter g 7 and that its generic asymptotics (g ^ 1) dif- 
fers from that of the finely tuned free fermion case g = 1 . 
Also, Appendices [A] and [B] provide details on the deriva- 
tion of the bosonized boundary theories, while Appen- 
dices [C] and [D] solve for the delocalized Majorana mode 
at the junction in the non-interacting and interacting he- 
lical/topological cases, respectively. 



II. TOPOLOGICAL 

SUPERCONDUCTOR-HELICAL LUTTINGER 

LIQUID JUNCTIONS 

The first junction we will analyze is that formed by 
a helical Luttinger liquid adjacent to a ID topological 
superconductor supporting a single localized Majorana 
mode at each end as illustrated in Fig. flTa). While ulti- 
mately we wish to understand the universal properties of 
the junction at low energies in the presence of (possibly 
strong) interactions in the Luttinger liquid, here we will 
begin by exploring the free fermion case which provides 
a useful point of reference. 



A. Scattering problem for free fermions 

Throughout this section we will assume that the su- 
perconductor is fully gapped save for the end-Majorana 
modes 71 and 72 shown in Fig. [IT a). We further as- 
sume that the superconductor is sufficiently long that 
these Majoranas overlap negligiblySEO anc i that the heli- 
cal wire couples only to 71 . The only relevant low-energy 
degrees of freedom are then 71 and those of the gapless 
helical wire, taken here to be non-interacting. At low 
energies it suffices to linearize the kinetic energy for the 
helical wire and incorporate the effects of the supercon- 
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FIG. 1: (a) Topological superconductor forming a junction 
with a helical Luttinger liquid. We assume that the helical 
Luttinger liquid couples only to the Majorana 71 at the inter- 
face, while its partner 72 remains a robust zero-energy mode. 
(b) Flow diagram for the junction as a function of the interac- 
tion parameter g for the Luttinger liquid. Provided g > 1/2, 
coupling to the Majorana 71 causes the Luttinger liquid to 
flow onto a fixed point where perfect Andreev reflection occurs 
at the junction. Here the zero-bias conductance for the junc- 
tion is quantized at 2e 2 /h. For Luttinger liquids with strong 
repulsive interactions such that g < 1/2, however, coupling 
to 71 is irrelevant. Here the system flows onto a fixed point 
where the Luttinger liquid undergoes perfect normal reflection 
at the junction, leading to a vanishing zero-bias conductance. 



ductor through local perturbations to the wire Hamilto- 
nian acting at the interface. (Formally, the latter pertur- 
bations can be derived by coupling the wire and super- 
conductor at the junction and then integrating out the 
gapped superconductor degrees of freedom.) Taking a 
semi-infinite wire at x > 0, we write the full Hamiltonian 
as H = Hq + 5H , where 5H represents the boundary 
terms and the kinetic energy Hq reads 



#o 



dx 



(-ivFipRdxi/iR + iv F ip L d x ip L ) . (1) 



Here vf is the Fermi velocity and ijr R , L creates right/left- 
moving excitations near the Fermi energy. It is conve- 
nient to rewrite Hq in terms of a single fermion field ip{ x ) 
defined over all x as follows: 



ip(x) 



i/)r{x), x > 
iPl(—x), x < 0. 



In terms of ip(x), Hq becomes simply 

Hq = dx (—ivFtp'd x ipj 



(2) 



(3) 



The boundary Hamiltonian 5H will encode various 
couplings at the junction, including tunneling between 



the Majorana mode 71 and the wire, local chemical po- 
tential modifications, and Cooper pairing induced locally 
in the wire by the superconductor. Thus the full Hamilto- 
nian will take the form of a Bogoliubov-de Gennes equa- 
tion. Before turning to a specific form of SH it is useful 
to discuss the problem in some generality. It is natural 
to expect that in the low-energy limit the effect of these 
couplings will be to impose certain boundary conditions 
on the wire, the precise nature of which will depend on 
details of the Hamiltonian. The set of possible bound- 
ary conditions corresponds, in renormalization group lan- 
guage, to different fixed points whose stability we would 
like to understand. For free fermions this issue can be 
exhaustively addressed using scattering theory. Here, all 
information regarding these putative fixed points as well 
as their stability is encoded in the S-matrix S(E), which 
at a given energy E relates the particle and hole states 
incident on the superconductor with the reflected states. 
Interestingly, using arguments similar to Refs. l48l49l one 
can show very generally that the S-matrix for an infinite 
wire can take only one of two forms in the E — > limit, 
corresponding to fixed points describing either perfect 
normal reflection, or perfect Andreev reflection at the 
junction. In the case of a finite wire of length L, the 
S-matrix has been derived in Refs. [27 50 51 , and one 
recovers the above answers in the appropriate limit, i.e. 
I = oo,£^0. 



1. General properties of the S-matrix 

The full Hamiltonian H including boundary couplings 
is diagonalized with quasiparticle operators carrying en- 
ergy E of the form 



T E = <f>Eli + / dxe 1 "f [P E (x)4,{x) + H E (x)^(x)}, 

J — OO 

where (f>E is the component of the wavefunction at the 
Majorana and Pe(x), He{x) respectively determine the 
particle- and hole-like amplitudes of the wavefunctions. 
Since Pe(— 00) and He (—00) describe left-moving elec- 
trons and holes incident on the superconductor, while 
Pe(oo) and He {00) capture the reflected right-moving 
states [see Eq. pj], the S-matrix is defined by 



Pb(oo) 

He(oo) 



= S(E) 



Pb(-oo) 

Ps(-oo) 



(5) 



Spp(E) S PH (E) 
Shp(E) Shh(E) 



Pe(-oo) 
He(-oo) 



To restrict the form of S(E), we first note that Te = 
r_ B , which follows from particle-hole symmetry exhib- 
ited by any Bogoliubov-de Gennes equation, combined 
with the fact that the helical wire has only a single 
fermionic species. This relation connects the particle 
and hole wavefunction amplitudes via Pl E {x) = He{x), 



which in turn implies that the S-matrix must obey 
S(E) = <T X S*(—E)<7 X , where a x is a Pauli matrix. Thus 
it suffices to determine (say) Spp(E) and Spn(E), since 
the other matrix elements follow from 



Shh(E) = S*pp(—E) 
Shp(E) = S PH (—E). 



(0) 



At zero energy Eqs. ^ allow one to restrict the S- 
matrix further still. Imposing Spp(0) = S^ H (0) and 
Spn{0) = Stf P (0) along with unitarity of the S-matrix as 
required by current conservation, one finds that S(E = 0) 
can take only two possible forms: 



S(0) 



5(0) 



o 











(perfect normal reflection) (7) 



(perfect Andreev reflection), 



( 8 ) 
for some phases a and ft. The diagonal form in Eq. 

n\ corresponds to the situation where particles incident 
on the superconductor undergo perfect normal reflection 
at the interface — an electron reflects as an electron with 
unit probability and similarly for holes. Conversely, the 
purely off-diagonal form in Eq. (l8| describes perfect An- 
dreev reflection at the junction; here electrons scatter 
perfectly into holes and vice versa, transmitting a Cooper 
pair into the superconductor in the process. In renormal- 
ization group terms, these limits correspond to two differ- 
ent fixed points at which the topological superconductor 
imposes either perfect normal reflecting or perfect An- 
dreev reflecting boundary conditions on the helical wire. 
Given our assumptions fixed points with intermediate 
boundary conditions are not possible (at least for free 
fermions) . 

Following Ref. [2^1 these two possible fixed points can 
be distinguished by the conductance across the junction. 
At a bias voltage V, the current transmitted into the 
superconductor is given by 



' = x 



dE\S PH (E)f 



(9) 



where \Sph(E)\ 2 is the probability that an incident elec- 
tron at energy E Andreev reflects into a hole at the 
junction, transmitting charge 2e into the superconduc- 
tor. The differential conductance G = 4h at T = is 
then 



G 



2e 2 



'Pfl 



eV)\\ 



(10) 



which in the zero-bias limit becomes 



G(V -> 0) = 



2r J 

h 



(perfect normal reflection) 
(perfect Andreev reflection). 



(11) 



While it is intuitively clear that perfect normal reflection 
ought to give rise to a vanishing zero-bias conductance, 
it is interesting that in the Andreev limit one necessarily 
obtains conductance quantization. We will return to this 
issue below when we obtain the S-matrix for a specific 
form of SH. 



2. Accessing the perfect normal and Andreev reflection fixed 
points with free fermions 

We would like to now understand the conditions re- 
quired for our non-interacting helical wire to flow onto 
each of the two possible fixed points identified above. 
With this objective in mind we will now consider the 
following boundary Hamiltonian: 



5H = 



dx 



t 



V2 
+ (iAipd x tp + h.c.) 



7i(V> — V-') + 2ui/j*i/j 

6(x). (12) 



Here, t allows electron tunneling between the Majorana 
7! and the helical wire at the junction, u is a local poten- 
tial that favors normal reflection, and the A term (which 
must involve a derivative by Fermi statistics) encodes 
processes wherein a Cooper pair hops between the helical 
wire and the superconductor. Additional couplings are in 
principle present but necessarily carry higher derivatives 
than those already displayed and can thus be safely ne- 
glected at low energies. Given this form of 5H, a straight- 
forward solution of the S-matrix yields, 



Spp 
Sph 



E[(u + i) 2 + A{i 2 + E 2 A)} 



it 2 + E[l + u 2 + A{t 2 
i(i 2 + 2E 2 A) 



E 2 A)} 



ii 2 +E\l 



? ,z 



A(i 2 + E 2 A)] 



(13) 



where the tildes denote quantities expressed in units of 
vf, e.g., u — u/vf- [Recall that the other two matrix 
elements follow from Eqs. dfjj.] For generic values of 
the couplings one clearly sees that the S-matrix becomes 
purely off-diagonal in the E -^ limit, corresponding 
to the onset of perfect Andreev reflection at the bound- 
ary and a quantized zero-bias conductance of 2e 2 /h for 
the junction. Avoiding this outcome requires fine-tuning, 
indicating that the Andreev fixed point is stable in the 
non-interacting case. 

Closer inspection of Eqs. ( 13 1 reveals that the flow 



towards perfect Andreev boundary conditions originates 
exclusively from the coupling to the Majorana 71 located 
at the junction. Indeed, in the E — > limit the S-matrix 
becomes purely diagonal upon fine-tuning t = 0, so that 
without coupling to the Majorana the system flows in- 
stead onto the (unstable) perfect normal reflection fixed 
point. (One might naively expect that the pairing term 
A alone would be sufficient to drive the system to the An- 
dreev fixed point, but this is not the case. In the helical 



wire, Pauli exclusion, which is responsible for the deriva- 
tive in this term, renders this an irrelevant perturbation 
at the perfect normal reflection fixed point. This can be 
understood by noting that when E — > 0, A drops out en- 
tirely from the S-matrix.) It follows that here the quan- 
tized 2e 2 /h conductance for the junction at the Andreev 
fixed point reflects the familiar zero-bias anomaly associ- 
ated with tunneling onto a Majorana mode 2 ^ 30 ! 32 ! 33 * 35 ] 
The scattering approach adopted here in fact follows 
closely the treatment of Ref. [29] who showed that reso- 
nantly coupling a non-interacting system to a Majorana 
mode generically induces perfect Andreev reflection as 
recovered here. One can intuitively understand this by 
observing that for any non-zero t, the Majorana zero- 
mode originally described by 71 gets absorbed into the 
helical wire, where it becomes a delocalized plane- wave 
state (see Appendix [C| . The Majorana character of this 
plane- wave must be preserved, however, which in turn 
guarantees perfect Andreev reflection. 

Next we explore the stability of the fixed points cap- 
tured here for free fermions when interactions are present. 
Interestingly, we will show that strong repulsive interac- 
tions modify the physics of the junction qualitatively. We 
will first treat the perfect normal reflection fixed point 
and then turn to the Andreev fixed point. 



B. Stability of the perfect normal reflection fixed 
point with interactions 

In the non-interacting limit accessing the perfect nor- 
mal reflection fixed point required fine-tuning to zero the 
coupling between the helical wire and the Majorana 71 
at the junction. Guided by this case, we will initially 
neglect the presence of ji and derive a fixed-point action 
describing perfect normal reflection for the Luttinger liq- 
uid. This is conveniently achieved using bosonization, 
where the right/left-moving fermionic modes in the wire 
are expressed in terms of dual bosonic fields <f>, 9 via 



junction: 






(14) 

(15) 



Physically, 9 relates to the fermion density n in the Lut- 
tinger liquid according to n — d x 9/Tt, while (f> and n 
are canonically conjugate variables. It will be useful be- 
low to observe that since one can write 9{x) — 9(0) = 
— 7rJ dx'n(x'), the fermion parity in the helical Lut- 
tinger liquid is given by Pll = cos[9(x = L) — 9(x = 0)], 
where x — L corresponds to the right endpoint of the 
wire. 

To ensure perfect normal reflection at the junction the 
fermionic fields are constrained to satisfy 



iP R {x = 0) = e la Tp L {x = 0) 



(16) 



for some unimportant phase a which we will simply set 
to zero. This in turn implies pinning of the field 9 at the 



0(x = 0) = mod 7T. 



(17) 



Similar constraints of course apply to the right end of 
the wire at x = L, so for concreteness we will henceforth 
set 9(x = L) = (except in Appendix [d] where a differ- 
ent convention is specified) . The fermion parity in the 
Luttinger liquid then reduces to 



P LL = cos[9(x = 0)]. 



(18) 



Thus the two pinning values in Eq. (17 1 correspond to 



the cases where the Luttinger liquid accommodates an 
even and odd number of electrons. 

Let us now arbitrarily select a particular pinning value 
for 9(x — 0). One can obtain an effective action for 
the remaining fluctuating phase field at the boundary, 
(f)(x = 0), in the following manner. First, the kinetic 
energy in Eq. (fTl) supplemented by interactions in the 
helical wire bosonizes to 



#o 



dx^[g{d x 4>f 

Z7T 



g-\d x 9f 



(19) 



where vp is the Fermi velocity and g is the Luttinger pa- 
rameter specifying the interaction strength (again, g = 1 
is the free-fermion limit while g < 1 and g > 1 re- 
spectively correspond to repulsive and attractive inter- 
actions). Obtaining the Euclidean action corresponding 
to Eq. ( 19 ) and integrating out all fields away from x = 
(see Appendix [A| for details) then leads to the following 
action, 



S 



normal 



dui 
2,1 / 2tt' 



Ml*l 



(20) 



where $ = <j>(x = 0). This action describes the perfect 
normal reflection fixed point, whose stability we can now 
assess. 

The most relevant perturbation to Sxiorm&\ originates 
from coupling to the neglected Majorana mode 71, 

ss t = ±=JdT 11 [Mx = ^-Mx = $)\ (21) 

which promotes Andreev processes at the junction. [The 
use of ipR as opposed to i\>l here is immaterial be- 
cause of the boundary condition of Eq. (16).] Bosoniz- 



ing 5 St requires some care, as the us ual procedure of 



naively replacing ipR ~ e in Eq. (21 1 leaves one with 



a non-Hermitian operator (among other technical prob- 
lems stemming from the Majorana operator). Indeed, 
since Eq. (21 1 contains 7! our bosonization must include 



the topological superconductor as well in order to ob- 
tain consistent results. We will now show how this can 
be done by considering a lattice model that includes the 
relevant low-energy operators for both subsystems. 

Because the only low-lying degrees of freedom in the 
superconductor are 71 and 72, we can distill this part 



down to a single fermionic lattice site taken to lie at po- 
sition 0: 



7i = 4 + c 
72 =i(4-Co). 



(22) 



Even though we formally model 71 and 72 as deriving 
from the same site, physically these operators are spa- 
tially well-separated. For our purposes the only conse- 
quence of this is that the Luttinger liquid, which we now 
define on a lattice indexed by sites j > 0, couples only 
to the linear combination c + cq on site 0. We there- 
fore consider the following tunneling Hamiltonian that 
hybridizes this site and the Luttinger Liquid, 



6H t = *7i(c{ - Ci) = t(4 + c ) (4 - c i)- 



(23) 



The full Hamiltonian is then 



H = H LL +H U (24) 

with the Luttinger Liquid Hamiltonian 

H LL = -J^>tc,- +1 + h.c.) + H int , (25) 

where H mt encodes density-density interactions. 

Next we implement a Jordan- Wigner transformation 
to write the Hamiltonian in terms of hard-core bosons bj 
via 



exp 






J 3i 



(26) 



with rij = b;bj = C-Cj. Furthermore, we rewrite the 
bosons at site in terms of Pauli spin matrices via 60 = 
((T x + icr v )/2, and those on the remaining sites j > in 
terms of a phase field fa: bj ~ ie l ^ j . Using exp (iirno) = 



a z the tunneling term H t in Eq. ( 23 1 then becomes 



SH t = ta x cos $, 



(27) 



where $ = fa. The Euclidean action corresponding to 
Eq. p7] is 



6 S t = 2t / dra x cos®, 



(28) 



which upon trivially rescaling t is the correctly bosonized 
form of the Majorana tunneling term in Eq. (21 1. 

For concreteness, it is useful to relate the wavefunc- 
tions for the superconductor in the fermionic and 'spin' 
languages. When t = the superconductor admits two 
degenerate ground states with well-defined but opposite 
fermion parity, |0) and |1), due to the Majorana zero- 
modes 71 and 72. These states are connected by the 
Majorana operators: 7i|0) oc |1). In the spin language 
|0) and |1) are eigenstates of a z , but the natural pair 
of degenerate ground states is formed by eigenstates |+) 



and |— ) of a x . The distinction between the two can be 
sharpened by considering a more realistic model for the 
superconductor, consisting of a Kitaev modePwith many 
sites intervening between 71 and 72. The ground states 
of given fermion parity, natural in the original fermionic 
representation, are then Schrodinger cat states for the 
spins, i.e. linear combinations of the two phase eigen- 
states, in which all spins point along the plus or minus 
x directionPS (Incidentally this is essentially why a ID 
spinless p-wave superconductor is a widely sought topo- 
logical phase of matter while the Ising spin chain is not, 
despite the fact that these models are superficially re- 
lated.) 

We are now in position to analyze the stability of the 
perfect normal reflection fixed point. The scaling dimen- 
sion of cos $ is l/(2g) while that of a x is zero at this fixed 
point, so under renormalization t flows according to 



dt 
d£ 



[1 - (29)"% 



(29) 



Equation ( 29 1 determines the renormalized coupling 



strength at a length scale I in terms of £ = ln(l/lo), with 
Iq a microscopic length of order the Fermi wavelength. 
Remarkably, for helical Luttinger liquids with g < 1/2 
tunneling onto 71 thus constitutes an irrelevant pertur- 
bation; perfect normal reflection is then stable despite 
the presence of a zero-energy Majorana mode to which 
the system can couple. We will explore the physical con- 
sequences of this result in the discussion. This coupling 
is relevant, however, when g > 1/2, indicating instabil- 
ity of the perfect normal reflection fixed point (consis- 
tent with our scattering analysis for free fermions with 
g = 1). Since this perturbation promotes Andreev reflec- 
tion at the junction it is natural to expect that at low en- 
ergies normal reflection then becomes entirely suppressed 
in favor of Andreev processes, just as we found for free 
fermions. We establish in the next subsection that this 
is indeed the case by examining the stability of the per- 
fect Andreev reflection fixed point when interactions are 
present. 



C. Stability of the perfect Andreev reflection fixed 
point with interactions 

One illuminating method for extracting the fixed point 
action describing perfect Andreev reflection at the junc- 
tion is to apply a duality transformation to the bosonized 
action S norma \ + SS t - In particular, this allows us to ex- 
tract the most relevant operator around this fixed point 
and analyze its stability. Our starting point is the parti- 
tion function for the perfect normal reflection fixed point 
expressed as a path integral; this is carefully derived in 
Appendix [B] and reads 



Z= T>$ 



r»=±l 



-2i/d-r<r x cos* 



(30) 



An important point here is that there is only a single 
sum over a x = ±1 (rather than one at each imaginary 
time slice). This arises from the fact that a x is a con- 
served quantity in the Hamiltonian and therefore has no 
imaginary time dynamics. Upon writing the cosine term 
in the Villain approximation the partition function then 
becomes 



Z 



x e 



™ E E « 

<7*=±1 {„( T )} e z 
* / <2T[$+(l+er x )7r/2-27™] 



S no „ 



(31) 



In Villainized form the t term can be decoupled with a 
Hubbard-Stratonovich field p, yielding 



Z 



x e 



V$>Vp ]T £ e~ s — 

<T a ==±l{n(r)}eZ 
/ d-r{ £ r +2ip[*+(l+£r x )7r/2-27m] j 



(32) 



At this point it is convenient to introduce the variable O 
which is canonically conjugate to $/7T by writing 



d T 
~2tt~' 



(33) 



Performing the sum over integers n(r) in Eq. (32) re 



stricts the field O to elements of 7rZ, while the a x sum 
simply imposes periodic boundary conditions on O along 
the imaginary time direction [i.e., 0(r — f3) — 0(r = 0) 
mod 2ir]. As usual the former restriction is difficult to 
handle so we impose the constraint 'softly' by adding a 
potential to the action which energetically favors integer 
values for Q/n; the partition function then reads 



X>$X>0e 



-Sn 



■ fdr 



-+ i -*d T e-t>cos(2e) 



(34) 



One can then integrate out $ to obtain the desired dual 
theory for 0: 



Z 



VOe 



-s dua 



(35) 



.5', 



dual 



t^l l 2 -W^<-<^>- <•><» 



where we dropped a term proportional to (<9 r 0) 2 since 
it is irrelevant compared to the first term in the action 
above. 

Since and $/7r are conjugate variables, the operator 
e 2 * e shifts $ by 2n. Thus the v term above represents 
an instanton operator which tunnels between adjacent 
minima of the <j x cos <& potential in Eq. ( 30 ) . [One might 
naively expect instanton operators that simultaneously 
change a x — > —<j x and $ — >• $ + it to be important, but 
these are forbidden since a x is a non-fluctuating classical 
degree of freedom.] Tunneling events imposed by cos(20) 



are qualitatively unimportant at 'large' t but must be 
retained otherwise. When t is relevant and the perfect 
normal reflection fixed point is unstable, the fixed point 
described by the dual theory with v — should therefore 
be stable and vice versa. This strongly suggests that 
Sduai with v — describes the perfect Andreev reflection 
fixed point. We will now confirm this by rederiving Sduai 
beginning from the fermionic theory. 

To access the perfect Andreev reflection fixed point, 
tunneling onto the Majorana 71 (which as we saw in Sec. 
II A 2 underlies the flow to this fixed point) must be in- 



corporated non-perturbatively. The coupling to 71 con- 
strains the fermionic fields at the junction such that at 
low energies 

$ t (x = 0) = e if> i/> L (x = 0), (37) 

which upon setting )3 = for simplicity pins the 
bosonized phase field <b to 



$ = (j)(x = 0) = or 7T. 



(38) 



(See Appendix [Cl for an explicit solution that derives 
these boundary conditions in the non-interacting limit.) 
The appearance of two possible pinning values for <f> 
can be understood from our analysis of the perfect nor- 
mal reflection fixed point above. There we showed that, 
by bosonizing the low-energy degrees of freedom for the 
Luttinger liquid and the topological superconductor, the 
coupling to 71 can be written as SH t ex ta x cos <&, where 
a x swaps between the two opposite-parity degenerate 
ground states |0) and |1) for the superconductor. Note 
that SH t commutes with the total fermion parity opera- 
tor 



Ptot =<r z cos9(x = 0), 



(39) 



with a z = 17172, as expected since the hopping preserves 
the global parity of the system. For g > 1/2 we showed 
that tunneling onto 71 constitutes a relevant perturbation 
at the perfect normal reflection fixed point. The system 
then flows at low energies onto a fixed point at which 
a x cos& is pinned to —1 (assuming t > 0). This allows 
two possibilities: 



a x = +1, cos* = -1 

a x = -1, cos* = +1. 



(40) 



corresponding to the two pinning values in Eq. ( 38 ) 



It is instructive to examine the ground states corre- 
sponding to the two sectors identified above. At the An- 
dreev fixed point these may be written as 



t) 
I) 






1) 



-l) 



|$=7T) 

8|* = 0), 



(41) 



where \a x = ±1) = |0) ± |1), and |$ = 0,tt) are the 
ground states of the bosonized Hamiltonian in Eq. (19 1 
with boundary conditions <p(x = 0) = 0,7r. Appropriate 
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linear combinations of the states in Eqs. (41 I yield ground 
states with well-defined parity Ptot, 



t> - U) 



(42) 



satisfying P to t|±) = d=|d=) . While the overall parity Ptot 
is a good quantum number, the individual parities of the 
superconductor and the Luttinger liquid are entangled in 
the two ground states |±). Explicitly, they can be written 
as 



(43) 



Here |Pll = ±1) = |$ = 7r) ± |$ = 0) have well defined 
parity in the Luttinger liquid. We note that the decom- 
position ( 43 1 is only exact in the limit of a semi- infinite 



|0) <g> |P LL = 


= -1) + |1) ® |P LL = 


= 1 


|0) <g> |P LL = 


= 1) + |1) ® |P LL = 


-1 



Luttinger liquid and when the Andreev fixed point has 
been reached. 

In fermionic langauge, the leading perturbation away 
from the Andreev fixed point is the potential term 



SS U — 2u dr 



ip R (x = 0)V>l(> = 0) + h.c. (44) 



which promotes normal reflection at the junction. Be- 
cause of the Andreev boundary condition ip R (x — 0) = 
ipt{ x = 0) the definition of 5S U above is somewhat sub- 
tle. When regularized at length scale e (for example via 
point-splitting) the operator ip R (x = 0)ijji,(x — 0) scales 
to linearly with e, so a divergent factor of 1/e must be 
absorbed in u to obtain a finite result. The bosonized 
form of this finite coupling is 



5S u ~4u J dr cos(29), 



(45) 



which is just the instanton operator in the dual action of 
Eq. (36) obtained earlier by complementary means. 



Since cos(20) has scaling dimension 2g, to leading or- 
der the coupling u flows according to 



du 
d£ 



(1 - 2g)u. 



(46) 



Provided g > 1/2 this perturbation is therefore irrelevant 
and the Andreev fixed point is stable. Note that this 
regime includes the non-interacting limit, g — 1, consis- 
tent with our scattering analysis above. For g — 1/2 the 
perturbation turns out to be exactly marginal: there ex- 
ists a line of RG fixed points which interpolate between 
normal and Andreev reflection. 

For a strongly repulsive wire with g < 1/2, however, 
the potential u is relevant. The Andreev fixed point is 
then unstable towards the perfect normal reflection fixed 
point analyzed previously (which we found is stable in 
this range of g). With g < 1/2 the backscattering term 
in Eq. (45 1 pins cos(20) to —1 at low energies (assuming 



u < for concreteness), yielding the two possible pin- 
ning values of 9{x — 0) identified at the perfect normal 
reflection fixed point in Eq. (17 1. Because of the bound- 
ary condition of fixed 8(x = L) these two pinning values, 
corresponding to different fermionic parities on the wire, 
are not degenerate, and in fact are split by an energy of 
order 1/L. 

It is interesting to ask about the fate of the Majo- 
rana zero-mode at the junction in this strongly repulsive 
regime. Although the tunneling t between 71 and the 
Luttinger liquid is irrelevant for g < 1/2, this term nev- 
ertheless has a quantitative effect on the zero-mode oper- 
ator since any finite t makes the commutator [H, 71] non- 
zero. Let 7 1 iew be the operator for the Majorana mode 
that arises when t / 0. In principle this operator can 
be determined by requiring that [H, 7^ ew ] = 0. We will 
alternatively deduce the asymptotic form of the probabil- 
ity density P{x) associated with 7"°™ using scaling. [It 
is tempting to view P(x) as deriving from the Majorana 
wavefunction corresponding to 7" ew — which is certainly 
legitimate in the free-fermion limit but rather subtle in 
the interacting case. In the regime of interest here with 
g < 1/2, 7" cw does not generally admit a single-body ex- 
pansion in terms of microscopic fermion operators. Ref- 
erence 52 discusses this issue and demonstrates that even 
for a strongly interacting system a Majorana 'wavefunc- 
tion' yielding a probability P{x) can be extracted from 
matrix elements of a fermion operator at position x with 
respect to opposite-parity ground states.] 



Dimensional analysis together with Eq. (29) lead to 
the following scaling ansatz for P(x), 



P(x) = ^p(x/L;tL 1 -^r 1 ^ 



(47) 



where L is the length of the Luttinger liquid. Note that 
V is a symmetric function of t that vanishes when t = 0. 
For g < 1/2 it suffices to treat t perturbatively since 
this coupling is irrelevant; at lowest nontrivial order one 
obtains 



P(*) ~ \ 



tL 



l-(29)" 



2 V{x/L). 



(48) 



In the perturbative regime we expect the Majorana wave- 
function to be normalizable in the L — > 00 limit (as op- 
posed to plane-wave- like as it is when t is relevant). Un- 
der this assumption the probability density P(x) must 
be independent of L in the thermodynamic limit which 
requires the asymptotic form 



V{x/L) 



(!)' 



-i/s 



(49) 



As an important self-consistency check, we note that for 
g = 1/2 (where t is marginal) Eq. (49) yields P(x) ~ 
1/x so that the Majorana wavefunction is only quasi- 
normalizable. For any g below 1/2, however, one obtains 
a normalizable probability distribution consistent with 
our ansatz. Thus we conclude that in the g < 1/2 regime 
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the irrelevant tunneling between 71 and the Luttinger 
liquid results in a Majorana zero-mode that bleeds into 
the wire but remains power-law localized to the interface. 
It would be very interesting to test these predictions in 
DMRG simulations by adapting the techniques of Ref. |55] 
to 'see' the Majorana wavefunction in this geometry in a 
numerical experiment. 

We also note as an aside that the critical value of 
g = 1/2 coincides with the critical Luttinger parameter 
below which pairing induced by proximity in the bulk of 
a helical Luttinger liquid is an irrelevant perturbation 53 . 
Thus helical wires with g < 1/2 not only resist Majorana 
modes imposed externally by a topological superconduc- 
tor as in the junction studied here, but also reject pairing 
that would make the wire itself topologically supercon- 
ducting. 



III. NON-TOPOLOGICAL 

SUPERCONDUCTOR-HELICAL LUTTINGER 

LIQUID JUNCTIONS 

We turn next to junctions formed by ordinary su- 
perconductors adjacent to helical Luttinger liquids, as 
sketched in Fig. |2ja). One may intuitively expect that 
here all traces of Majorana physics captured in the pre- 
vious section are simply absent, since the superconduc- 
tor no longer supports protected Majorana zero-modes. 
While this is indeed true in the non-interacting limit, we 
will demonstrate that Majorana modes can generically 
appear in the case of interacting helical Luttinger liq- 
uids. To establish this counterintuitive result, we first 
note that in our analysis of generic properties of the S- 
matrix for free fermions in Sec. |II A 1| the topological 
nature of the superconductor we were considering there 
played an irrelevant role. The same analysis [with the 
sole modification of dropping the </>e7i term in Eq. Q] 
applied to present case leads to identical conclusions — in 
the E — > limit the S-matrix must again be either purely 
diagonal or purely off-diagonal. Thus even in the case of 
a junction formed with an ordinary superconductor, two 
physically allowed fixed points for free fermions remain, 
corresponding to perfect normal reflection and perfect 
Andreev reflection at the interface. We will adopt a sim- 
ilar program to that followed in the previous section, be- 
ginning by understanding how to access these two fixed 
points for free fermions, and then addressing the stability 
of these fixed points when interactions are present using 
bosonization. 



(a) 



Ordinary SC 



Helical LL 



(b) 




Ordinary SC fflllS Helical LL 

sir 

/ 



(c) 



Perfect Andreev Reflection 



u w v w;v v u u A A A A 
\i \t \t w'v u v u A A A A 



g=0 g=1 g=2 

Perfect Normal Reflection 

FIG. 2: (a) Ordinary superconductor forming a junction with 
a helical Luttinger liquid. If the superconductor is fully 
gapped throughout, then in the non-interacting limit the he- 
lical wire necessarily flows to a fixed point at which perfect 
normal reflection occurs at the junction. Accessing the per- 
fect Andreev reflection fixed point in the non-interacting limit 
requires the presence of a zero-energy Andreev bound state 
/ = (71 + *72)/2 at the junction as shown schematically in 
(b), with fine-tuning such that the wire couples only to (say) 
71. (c) Flow diagram for the junction as a function of the 
interaction parameter g for the Luttinger liquid. For helical 
Luttinger liquids with g < 2, the perfect normal reflection 
fixed point is stable. When g > 2, however, this fixed point 
is unstable towards the Andreev fixed point. Remarkably, 
here the Andreev bound state required to achieve perfect An- 
dreev reflection will be generated dynamically, and with no 
fine-tuning required. 



full Hamiltonian for the problem as H = Hq + SH, where 
Hq is the kinetic energy for the wire defined in Eq. (p| [or 
equivalently Eq. ([II] and SH contains the terms acting at 
the boundary. Retaining only the leading potential and 
pairing terms at the junction, one has 



Perfect normal and Andreev reflection fixed 
points for free fermions 



SH 



dx [2uip f ip + (iAiPd x ip + h.c.)] 8(x), (50) 



Let us start by making the (physically reasonable) 
assumption that the ordinary superconductor is fully 
gapped throughout, so that the helical wire hosts the 
only low-energy degrees of freedom. We again write the 



which simply corresponds to Eq. ( 12 ) considered previ 



ously without the Majorana term. The S-matrix for our 



Hamiltonian can therefore be read off from Eqs. (13 1 by 
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simply setting t = 0; this yields 



Spp 

Sph 



(u + i) 2 + {EA) 2 
l + u 2 + (EA) 2 

2iEA 
+ u 2 + {EA) 21 



(51) 



where again the tildes denote quantities normalized by 
v F . 

Notice that when E — > 0, A drops out and the S-matrix 
takes the form of Eq. (I7| corresponding to perfect nor- 
mal reflection. No fine-tuning of parameters in SH can 
alter this conclusion, nor can additional local perturba- 
tions involving only ip, ip' (which necessarily carry ad- 
ditional derivatives than those already included above, 
and are thus unimportant at low energies). Clearly then 
the perfect normal reflection fixed point is stable in the 
non-interacting limit. 

How, then, can one access the physically allowed con- 
dition of perfect Andreev reflection? For free fermions, 
accessing this fixed point requires the presence of an addi- 
tional localized Andreev bound state in the superconduc- 
tor, with corresponding operator /, to which the helical 
wire can couple; see Fig. (2jb) for a schematic illustra- 
tion. (Our conclusions for the general properties of the 
S-matrix discussed in Sec. |II A 1| again remain unaltered 
by the addition of this mode.) To see how this tran- 
spires, let us write / in terms of Majorana operators 71,2 
via / = (71 + «72)/2 and consider the following interface 
Hamiltonian, 



SH 



dx 



+ i^lil2 



t t' 

— 7 1 (v> t - VO + i—ppzw + vo 



5(x) 

(52) 



The t, t! terms encode the most general bilinear couplings 
between ip,ip' and /, p (after shifting these operators 
by overall phases to make the couplings real), while S 
sets the energy for the Andreev bound state, which will 
generically be non-zero. We have dropped the u and A 
terms considered above for simplicity since these consti- 
tute qualitatively unimportant perturbations here. With 
this form of SH the S-matrix elements are now given by 



Spp 
Sph 



(it 1 - iS) 2 



E 2 



{tt') 2 + S 2 -i(t 2 + t' 2 )E-E 2 
-i{i 2 - t' 2 )E 
{it') 2 + 5 2 - i{i 2 + i' 2 )E - E 2 ' 



(53) 



with t = t/vp. etc. in our usual notation. As before, one 
finds that in the limit E — > the off-diagonal elements 
still generally vanish, reflecting stability of the perfect 
normal reflection fixed point even in the presence of an 
additional Andreev bound state at the junction. One 
can, however, now fine-tune couplings in SH to obtain 
the desired perfect Andreev reflection at zero energy by 



taking 5 = and either t = 
ness let us choose t' = S = 
components become simply 

Spp = 
Sph = 



- or t! = 0. For concrete- 
0, upon which the S-matrix 

E 

it 2 + E 

it 2 



it 2 + E 



(54) 



Taking E — > in this fine-tuned limit produces an S- 
matrix of the form in Eq. (fSJ) , which indeed corresponds 
to perfect Andreev reflection. 

More physically, in the non-interacting limit perfect 
Andreev reflection in the E — > limit can emerge in this 
type of junction only when a pair of zero-energy Majo- 
rana modes appears at the boundary, with the helical 
wire coupling to only one of these. Thus precisely as in 
the topological superconductor junction, the absorption 
of the Majorana mode 71 into the helical wire underlies 
the onset of perfect Andreev reflection. One crucial dif- 
ference of course is that this Majorana's partner, 72, is 
spatially separated in the topological case, which allows 
the Andreev fixed point to be stable for free fermions. In 
the present context the coupling t' to 72 in Eq. ( 52 1 is 



generically non-zero, and constitutes a relevant perturba- 
tion which drives the system back to the perfect normal 
reflection fixed point. 

A second, and intimately related, difference relative to 
the topological case is that here 71 and 72 generally com- 
bine to form an Andreev bound state with finite energy S. 
Interestingly, provided t' = perfect Andreev reflection 
nevertheless survives when incident electrons in the heli- 
cal wire are resonant with this bound state. This can be 
understood from Eqs. ( 53 1 by noting that when t' = the 
S-matrix becomes purely off-diagonal at energy E = S. 
Similar behavior has been captured previously by Law et 
al¥** in the context of tunneling into discrete Majorana 
edge modes in a two-dimensional p + ip superconductor. 
However, we believe this behavior is special to the free 
Fermi case g = 1 and does not survive in the interacting 
context. 

Before turning to the interacting case it will prove ben- 
eficial to discuss the t' and S perturbations from a slightly 
different perspective. At the perfect Andreev reflection 
fixed point, the coupling t — which must be treated non- 
perturbatively — strongly constrains the behavior of ji at 
low energies. More precisely, at energies E -C t the dy- 
namics of 71 are 'slaved' to those of \jj and ^ such that 
the system avoids paying a large energy cost from the t 
term in Eq. ( 52 1 . One can show this explicitly by diago- 
nalizing the Hamiltonian H = Hq + SH with t! = S = 0, 



H = dx 



-ivf 



^d x ip 



- + —^-fi(^ - ^)S(x) 



(55) 



for a finite-size wire and then expanding ip and 71 in 
terms of the resulting low-energy modes. This calculation 
is sketched in Appendix O and yields the familiar perfect 
Andreev reflection boundary condition, ip(x = + ) = 
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i[>r(% = 0) = ip^{x = ) = ip L (x = 0), along with the 
following constraint on 71 at low energies: 

7i - *^=^ t (° + ) ~ ^(° + ) - ^(0-) + ^(0-)]. (56) 

One can understand this result more intuitively in the 
following manner. Let us decompose ip(x) in terms of 
two Majorana modes by writing 



tp(x) = 



7 r (x) +in/i(x) 



(57) 



in this new basis the Hamiltonian reads 
H = —i I dx 



(58) 

While j r (x) is unaffected by 71, the Majorana ji(x) ex- 
hibits a sign change at x — due to the coupling t. Ulti- 
mately this sign change gives rise to the perfect Andreev 
reflection boundary conditions. Since ji(x) behaves sin- 
gularly at the interface it is helpful to introduce a new 
Majorana operator ji(x) which is well-behaved every- 
where by defining 



ji(x) = sgn(a;)7 l (a;). 



(59) 



Differentiating the Euclidean action corresponding to Eq. 
( 58 1 with respect to 7^ yields the following equation of 



motion: 



= (d T — ivFd x )"ii + iv2t^/i5(x). 



(60) 



Rewriting this using Eq. ( 59 1 yields a second term pro- 



portional to S(x) due to the singular nature of the trans- 
formation at x = 0: 

= sgrx(x)(d T ji- WFd x *n)+i\/2t 71 —^7, J S(x). 

.( 61 ) 
At energies E <C t, the fields must conspire to eliminate 

the boundary terms above so that t disappears entirely 

from the equation of motion. This indeed occurs if we 

pin 71 = t " F 7i , which agrees with Eq. (56 1 derived by 

completely different means. 



It follows from Eqs. (52 1 and (56 1 that the two leading 



perturbations about the perfect Andreev reflection fixed 
point, given by turning on if and S, project to the same 
operator at low energies. 



B. Stability of the perfect normal and Andreev 
reflection fixed point with interactions 

Having understood how to access the physically al- 
lowed perfect normal and Andreev reflection fixed points 
for free fermions, we proceed now to assess the stabil- 
ity of each when interactions are present, beginning with 



the former. Just as in Sec. |IIB| bosonizing the Hamil- 
tonian and imposing perfect normal reflection boundary 
conditions pins 9(x = 0) to either or ir. Equation (20 1, 
repeated here for clarity, 



5. 



normal 



dui 

2n I 2^ 



WI4I 



(62) 



again describes the fixed point action for the fluctuating 
field $ = (f>(x = 0). Due to the generic absence of Ma- 
jorana modes at the junction, however, Eq. (21) no long 



provides the leading perturbation away from this fixed 
point. Instead, the most relevant perturbation which in- 
duces Andreev reflection at the junction corresponds to 
processes in which a Cooper pair hops from the helical 
Luttinger liquid into the ordinary superconductor: 



6Sa 



dr[iipUx = 0)ipl(x = 0) + h.c 



2 A / drsin(2$). 



(63) 



Equation ( 63 1 can be obtained via point-splitting in the 



unfolded chiral fcrmion theory: exp(2i<l>) is the leading 
operator appearing in the operator product expansion of 
two ^ fields. Since sin(2$) has scaling dimension 2/g, 
the coupling A flows according to 



rfA 
~dl 



(l-2<?- 1 )A 



(64) 



and is therefore irrelevant for g < 2. Thus the perfect 
normal reflection fixed point is stable not only for free 
fermions, but also for helical Luttinger liquids with arbi- 
trarily strong repulsive interactions, or attractive inter- 
actions below a critical strength. When attractive inter- 
actions exceed this critical strength, resulting in g > 2, 
A then constitutes a relevant perturbation which drives 
the system away from this fixed point. Physically, super- 
conductivity is induced at the endpoint and spreads to 
the rest of the wire, resulting in a 'topological' supercon- 
ductor, but without exponentially localized end states. 
Rather, because the superconductivity was seeded at only 
one endpoint, the Majorana modes both effectively live 
at that endpoint, with a splitting that turns out to be 
power law vanishing in the length of the wire. We will 
now argue more precisely that the helical Luttinger liq- 
uid flows to perfect Andreev reflection, and explain the 
physical consequences. 

Duality once again provides an effective tool for iden- 
tifying the fate of the system when g > 2. Following 



the steps outlined in Eqs. (30 1 through (36 1 to dualize 
^normal + <5Sa > one obtains the dual action 



°dual 



did 

2tt 2ng ' 



lei 



dr cos 0. 



(65) 



Here $ and <d/ir are once again conjugate variables, so 
that cos represents an instanton operator which tunnels 
between adjacent minima of the sin(2$) potential in Eq. 
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( 63 1 . Thus v' is expected to be irrelevant when A is rel 

with v' 



evant (and vice versa) , suggesting that S' dusd 



I) 



describes the physically allowed perfect Andreev reflec- 
tion fixed point that we identified in the free fermion case. 
Furthermore, since cos G represents the fermion parity in 
the helical Luttinger liquid, the dual action above also 
suggests the following: 1) at the fixed point with v' = 
there is a degeneracy between states with even and odd 
fermion number, and 2) the leading perturbations away 
from this fixed point split this degeneracy. We will now 
put these statements on firmer footing by explicitly con- 
structing the fixed point action and leading perturbations 
beginning from the fermionic theory. 

As in Sec. |IIC| the perfect Andreev reflection bound- 
ary condition pins <p(x — 0) = or 7r, and the fixed point 
action for O = 9(x — 0) derived by bosonizing the Hamil- 
tonian and integrating out the fields away from x — is 
given by 



5'a 



ndrccv 



du) 

2tt 1-ng ' 



ier 



(66) 



which indeed recovers S dual in the v 



limit. One 

might naively guess that the most relevant perturbation 
that induces normal reflection is the u term in Eq. ( 44 1 , 



but this is incorrect. To properly capture the physics it is 
crucial to recall how we accessed the Andreev fixed point 
in our scattering analysis for free fermions above. Doing 
so required the presence of two zero-energy Majorana 
modes at the junction, -fi and 72 in Fig.pfb), along with 
fine-tuning such that the helical wire absorbed 71 (say) 
but decoupled completely from 72. Without interactions 
perfect Andreev reflection at zero energy emerges only 
in this fine-tuned limit. One can then see that the lead- 
ing perturbations correspond to the 5 term in Eq. ( 52 ) 
which lifts the energy of 7^2 and the coupling between 
the Luttinger liquid and the Majorana mode 72, 



5S~ n = 



l^lll2 



V2 



dT 72 [tp R (x = 0) f + xj) R {x = 0)]. (67) 



Both perturbations were tacitly neglected by asserting 
that the system resided at the perfect Andreev reflection 
fixed point. We have explicitly verified in Sec. |III A| 
that both of these project to the same operator in the 
low energy limit, which corresponds to the one leading 
operator cos O that appears upon performing the duality 
transformation. 



Deducing the flow of v' as discussed in Sec. II B 
finds 



one 



dv' 

~dl 



[l-g/2]v'. 



(68) 



For g < 2 this coupling is relevant and drives the system 
back to the perfect normal reflection fixed point, which 
is indeed stable in this regime. More interestingly, for 



a helical Luttinger liquid with strong attractive interac- 
tions such that g > 2, coupling to this second Majo- 
rana mode represents an irrelevant perturbation, imply- 
ing stability of the perfect Andreev reflection fixed point. 
In other words the fine tuning required to access perfect 
Andreev reflection for free fermions is no longer necessary 
with strong attractive interactions. Figure p[c) summa- 
rizes the renormalization group flows found in this section 
for the ordinary superconductor/helical Luttinger liquid 
junction. Similarly to the case of the helical/topological 
junction, when g > 2 we expect one of the two dynami- 
cally generated Majorana modes to delocalize completely 
into the Luttinger liquid, and the other one to be power- 
law localized at the junction. 



IV. NON-TOPOLOGICAL 

SUPERCONDUCTOR— SPINFUL LUTTINGER 

LIQUID JUNCTIONS 

We now consider a junction between a non-topological 
superconductor and a semi-infinite spinful Luttinger liq- 
uid [see Fig. km a)]. In previous work on this problem^ 54 * 55 !, 
conductance and local density of states were calculated 
at the Andreev fixed point of such a junction. Here we 
emphasize that for weakly repulsive interactions such an 
Andreev fixed point is ultimately unstable, and the sys- 
tem generically crosses over to the normal reflecting fixed 
point. This instability of the Andreev fixed point is in 
fact crucial for establishing sharp transport signatures of 
Majorana modes in the topological case considered in the 
next section. 

We begin by bosonizing the problem. The two spin 
channels a =f, 4 are bosonized according to 



*pR/La 



= e (i<t> a ±ie a ) _ 



Defining the charge and spin fields 



(69) 



<t>p/<j = ("At ±<A|.), 

e p/a = (0 t ± e x ) /V2, (70) 

we then consider the general quadratic bulk action 



<Sbulk = 



2tt 



dx / dr 



-d x 6pd T 6 p + 

r 



-d x (j)ad x e 



2tt 



(dx^aY 



9T 

2tt 



(d x e p y 



9a 

2ir 



(d x o a y 



(71) 



Here g p / a are the Luttinger parameters characterizing in- 
teractions in the charge/spin sectors, and v the velocity. 
To facilitate a Luttinger liquid analysis, we first (implic- 
itly) assume a single velocity v and also g„ — 1, which 
corresponds to a spin-SU(2) invariant system. We then 
discuss how the analysis is modified in the experimen- 
tally relevant situation where SU(2) breaking terms are 
present. 
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(a) 



Ordinary SC 



Spinful LL 



(b) 

Perfect Andreev Reflection 
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FIG. 3: (a) Trivial superconductor forming a junction with a 
spinful Luttinger liquid, (b) Flow diagram for the junction as 
a function of the charge-sector interaction parameter g p for 
the Luttinger liquid, in the limit where the spin-sector inter- 
action parameter is g a = 1. Note that the free-fermion limit 
is very special — here there are marginal boundary couplings 
which lead to a non-universal zero-bias conductance ranging 
anywhere from to 4e 2 //i depending on parameters. For any 
repulsive interaction strength (g p < 1), however, the junction 
flows to the perfect normal reflection fixed point where the 
zero-bias conductance vanishes. 



There are two natural fixed points: normal reflecting, 
with 6 pi 9 a pinned at x = 0, and Andreev reflecting, with 
4> p and (/> CT pinned. The boundary action for the normal 
fixed point is obtained by integrating out all fields except 
for 4>p/ a {x — 0), resulting in 



s 



normal 



duj \uj\ 
2tt 2tt 


"Cflp + ftO (1*^ + 1*^) 


+ (3p-.9<x)$tM < M- w ) 



(72) 



where Q-\/\, = 4>f/i(x — 0). The leading perturbation to 
this boundary action is the Cooper-pair tunneling term 
cos(v / 2 < i ) p), which has dimension g^ 1 and is hence rele- 
vant for g p > 1. In this range of g pi the resulting Andreev 
fixed point t o w hich the system then flows is simply the 
dual of Eq. ([72J)P 



5' 



Andreev 



dui |w| 
2tt 2tt 



i% 1 +9^) 



(l©tl 2 + l©4| 2 ) 



+ (. 9p 1 -. g - 1 )e t Me i (- w ) 



(73) 



with Sf/4. = 9f/±(x = 0). The leading perturbation 
is now normal backscattering, described by cos(\ /f 2Q p ), 
with scaling dimension g p . 



The key point is that in the physically relevant regime 
of weak repulsive interactions, g p < 1, this normal 
backscattering term is relevant — i.e., the Andreev fixed 
point is unstable. Thus a junction between a spinful Lut- 
tinger liquid and an ordinary non-topological supercon- 
ductor is described by a stable normal reflecting fixed 
point. This stands in sharp contrast to the free fermion 
situation g p = 1, where these operators are exactly 
marginal and allow any value of the zero-bias conduc- 
tance between and 4e 2 //i, as can be seen by explic- 
itly solving the free fermion scattering problem. (This is 
actually a special case of the free fermion solution for 
the topological/spinful junction presented in the next 
section; indeed, upon setting the coupling to the Ma- 
jorana to zero one obtains the non-universal zero-bias 
conductance quoted here.) Arbitrarily weak repulsive in- 
teractions drive a crossover to the normal reflecting fixed 
point, where the conductance goes to zero. Figure p[b) 
summarizes the renormalization group flows for this case. 

Of course, in any proposed physical realization of a Ma- 
jorana wire, spia.SU {2) symmetry is broken by the Zee- 
man and spin orbit couplings (but in a simple Rashba 
model with density-density interactions and spin orbit 
coupling g a = 1 because of a hidden SU(2) symmetry^. 
At the level of free fermions, the only modifications these 
force on the low energy theory is differing Fermi momenta 
and velocities for the two species. The former is simply a 
restriction on the types of operators that can appear as 
perturbations (i.e. some operators might not be allowed 
because of a kp mismatch), while the latter appears di- 
rectly in the low energy action. To include interactions 
we must bosonize, and it is easiest to do so in a basis 
which diagonalizes the velocity; that is, we separately 
bosonize the two spin modes with the differing veloci- 
ties. The kinetic term then becomes an arbitrary 2 by 
2 symmetric matrix which generalizes the two Luttinger 
parameters g p and g a . In principle we can integrate out 
the bulk to obtain a boundary theory and analyze the rel- 
evance of Cooper pair tunneling as above. Although the 
resulting phase diagram depends in a complicated way on 
the three kinetic term coefficients and the two different 
velocities, it is still the case that the free fermion point is 
exactly marginal for any choice of velocities. Hence there 
is a robust region of interaction parameter space near the 
free fermion fixed point - roughly speaking the set of re- 
pulsive interactions - where we are driven to the perfect 
normal reflecting fixed point. Thus, even in the case of 
sp'mSU(2) symmetry breaking we generically expect no 
zero bias peak in the non-topological superconductor / 
spinful Luttinger liquid junction. 



V. TOPOLOGICAL 

SUPERCONDUCTOR— SPINFUL LUTTINGER 

LIQUID JUNCTIONS 

Finally, we analyze the junction sketched in Fig. 4(a) 
between a topological superconductor and a semi-infinite 
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spinful Luttinger liquid. As in the helical case we will 
first attack the free-fermion limit and then treat the in- 
teracting case using bosonization. 



A. Scattering problem for free fermions 



Here we have used time translation invariance to restrict 
to a solution with energy E. To handle the derivatives 
and delta functions at X = 0, we introduce the notation 
7i(0 ± ) = 7i(±e). At E = 0, the first equation of motion 
then gives 70 (0 + ) = — 70 (0~), and thus the remaining 
ones become 



We start with the free fermion scattering problem and 
compute the zero bias tunneling conductance when one 
allows for rather general boundary terms in the Hamil- 
tonian. Remarkably, we will show that this is robustly 
quantized at 2e 2 /h independent of parameters, so long 
as coupling to the Majorana zero-mode at the junction 
remains finite. Consider first the case when the spinful 
Luttinger liquid exhibits SU(2) spin rotation symmetry 
in the bulk. Coupling to the Majorana mode 7 at the 
junction can then always be gauge transformed to the 
fornP3 



8H X = iX dx 7 Oi + V4) 8{x), 



(74) 



where i\)\£ can be obtained from ipt,i v ^ a a rotation and 
a phase (we have replaced our previous label 71 with 7 
for future notational clarity). Neglecting terms involv- 
ing derivatives, the most general quadratic form of the 
boundary Hamiltonian at x — reads 

SH = J dxS(x)[i\j(ipi+iljl) 

+ {A 12 lpllp 2 + Vl2lp\lp2 + Ull^iV'l + V 2 2^l^2 + h.C.)] 

(75) 

Because A12 and v±2 arc complex while vn and ^22 are 
real, we have a total of 6 real parameters. The conduc- 
tance follows upon computing the S-matrix as a function 
of these 6 parameters and then summing up the proba- 
bilities of the Andreev processes. 

Let us now perform the calculation in detail. We first 
decompose ip\,2 in terms of Majorana fermion field oper- 
ators 7^ as follows: 



i>\ = 7o + «72 
1P2 = 7i + *73 



and rewrite Eq. ( 75 I as 



SH = 



dx 



A 7 7o(0) 



2 ^ 



3^7^(0)7^(0) 



(76) 



(77) 



Here g^ v is a real anti-symmetric matrix, encoding the 
6 free parameters introducted earlier. Let Sq be the free 
fermion bulk action of the spinful wire and SS be the 
boundary action corresponding to Eq. ( 77 1 . The equa- 
tions of motion that follow from varying So + SS are: 



= 



S{S + SS) 

Si 
S(S + SS) 

Sin 
iXjSfj.oSix) - 



= ^A7 (0) - E 1 



-idxln + i} j g^i„(0)S(x) 



El, 



(78) 



7i(0 + )=7»(0-)+^9«7i(0) (79) 

for j = 1,2,3. These can be summarized as 

7„(0 + )=E A V7,(0~) (80) 

where Man = -1, Mm = M,n = 0, and Mu = A~ 1 B, 



J oo 
with Aij — 5ij 



■9ij 






■9ij 



In fact, the 3x3 ma- 



trix Mij is just a generic element of 50(3) . Equation ( 80 1 
gives the scattering matrix in terms of Majorana fields. 
Changing basis to complex Fermi fields and calculating 
the four Andreev transmission probabilities tpj — > tpj , we 
find that they generically add up to 1; consequently the 
zero-bias conductance is indeed quantized at 2e 2 /h. 

What about the physically relevant SU(2) non- 
invariant case? For gapless free fermions, the only way to 
break ST/ (2) invariance is to introduce different velocities 
Vi 2 so that 



Ho 



dx 



E 

a=l,2 



\-iVa (i> ] Ra d x %l) Ra - i%l^ La d x il> 



(81) 
Now we cannot rotate the fields in such a way that only 
-01 couples to the Majorana mode 7, and are forced to 
retain both Ao77o(0) and Ai77i(0). Proceeding as in the 
St7(2)-invariant case we again derive Eq. (80 1, where now 
M = A~ l B. The matrices on the right side are given by 
A = V + G + L,B = V - G - L, with 4x4 matri- 
ces V, G, L defined as follows: V — diag(i>i,V2,«i,'U2), 
L^ v = —^Xjj,X v (where A M = for (jl > 1), and G M „ a 
general anti-symmetric ^-independent matrix containing 
the quadratic boundary couplings. 

In the case of differing velocities, the matrix M defined 
above is no longer equivalent to the scattering matrix, 
and in particular is not unitary. Indeed, the unitarity of 
the S'-matrix follows from the conservation of probability 
current, which is proportional to V. The correct S matrix 
is then given by 



S 



Vvm^= 

Vv 



Rewriting this equation as 

\W W 
and noting that 



1 B ^= 

v Vv 



1 A^- = 5, v + ^-{G + L) ^- 

V Vv M Vv VV 



(82) 



(83) 



(84) 
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FIG. 4: (a) Topological superconductor forming a junction 
with a spinful Luttinger liquid, (b) Flow diagram for the junc- 
tion as a function of the charge-sector interaction parameter 
g p for the Luttinger liquid when the spin-sector interaction 
parameter is g a = 1. For g p < 1/3 the system flows to a per- 
fect normal reflection fixed point characterized by a vanishing 
zero-bias conductance. When g p > 1/3, however, the junction 
flows to a novel fixed point corresponding to perfect Andreev 
reflection for one species and perfect normal reflection for the 
other, yielding a quantized 2e 2 jh conductance. 



and similarly for B, we can now perform the rotation 
along /i, v = 0, 1 to eliminate Ai. Taking the limit E — > 
we recover the same form of the S'-matrix as in the SU(2) 
non-invariant case, and the same quantized conductance 
G = 2e 2 /h. 

As remarked in the preceding section this result breaks 
down only when the coupling to the Majorana is fine- 
tuned to zero. In this special limit one obtains a 
free fermion non-universal zero-bias conductance rang- 
ing from to 4e 2 //i [in contrast to the helical case, Pauli 
blocking is absent here so that the A 12 term in Eq. (75) 



can efficiently transmit Cooper pairs into the supercon- 
ductor]. From this perspective it is somewhat curious 
that when the Majorana coupling is restored, A 12 is un- 
able to enhance the zero-bias conductance beyond 2e 2 /h. 
Evidently one of the channels [?/>i in the S'L r (2)-invariant 
limit] hybridizes with the Majorana mode 7, leading to 
Andreev boundary conditions but also blocking transport 
of ipiip2 Cooper pairs. We should emphasize that this re- 
sult is specific to having only two conducting channels. 



B. Phase diagram with interactions 

We now analyze the interacting case. We first bosonize 
as in the previous section, obtaining the bulk action in 
Eq. (71 1. Though we write expressions with general g a , 

i.e.. unbro- 



ken SU(2) spin rotation symmetry in the bulk of the 
Luttinger liquid. Then we proceed to discuss the physi- 
cally relevant 5/7(2) non-invariant case. As for the heli- 
cal/topological junction, we now need to include an ad- 
ditional effective spin- 1/2 degree of freedom for the Ma- 
jorana modes at the ends of the topological superconduc- 
tor. For simplicity we assume that the Majorana mode 
at the junction only couples to (say) the spin-up electron. 
Tunneling onto the boundary Majorana mode bosonizes 
to 



6S\ = A / dra x cos $ 



t- 



(85) 



[As before we define 3> T /4, — <f>^n{x — 0) and O^/i — 
9^n{x — 0)]. Having already dealt with the subtleties of 
the spin-1/2 degree of freedom in Sec. Ill] we now arbi- 
trarily fix the eigenvalue of a x to +1 and drop a x from 
the subsequent analysis. 

We begin with the normal reflecting fixed point real- 
ized at A = 0, where Q = 0. To determine its stabil- 
ity we have to co mpu te the scaling dimension of cos$-|-. 



which using Eq. (|72| is (g- 1 + g o r 1 )/4. With g a = 1, 



this term is relevant for g p > 1/3. In this case, the 
system flows to a novel Andreev/normal (A <g) N) fixed 
point where $^ and 0^ are simultaneously pinned. Here 
the spin-up electrons exhibits perfect Andreev reflection 
while the spin-down electrons undergo perfect normal 
reflection. Deriving the boundary field theory action 
for this new fixed point is somewhat subtle. We again 
start with the bulk action (71 1, but this time we inte- 
grate out </>^ and 0j,, taking care to respect the boundary 
conditions <f>f(x = 0) = 0^(x = 0) = (handling these 
boundary conditions incorrectly results in a spurious to- 
tal derivative term for the bulk that non-trivially modi- 
fies the boundary field theory, yielding the wrong scaling 
exponents). This yields 



2tt 



dxdr 



, I A d A) 2 + -1^ _i (W 2 



+ - / dxdr ^-^ ) {d T e^d x 4>i - d^drh) . 

W \9p + 9<yJ 

(86) 

Note that the second (Berry phase) term is a total deriva- 
tive, and can be integrated to give 



'Berry 



= — dr 



9 P - 9o 

9p + 9a 



e t 9 r $. 



(87) 



Integrating out 4>i (x) and (Mx) for x > from ( 86 ) and 
combining the result with ( 87 ) finally yields the following 
fixed-point action 



S 



A®N 



doj \ui\ 
2?r 2tt 



(ffZ + 'A 1 ) 



W 2 



{9p+9a) 



ie t | 2 



in the stability analysis we assume g a = 1, 



9p~ 9a 
9p + 9a 



dr#ia r e 1 
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From Eq. ( 88 1 we can determine all of the scaling di- 
mensions. In particular, ifi^d x ip^ ~ e 2 f has scaling di- 
mension 4/(g~ 1 +g~ 1 ), and is thus relevant precisely for 
g p < 1/3. In that case the system flows back to the nor- 
mal reflecting fixed point, whereas for g p > 1/3 it is the 
leading irrelevant operator around the new A <g) N fixed 
point. These renormalization group flows are summa- 
rized in Fig. Pub). 

We now perform a more general analysis, involving all 
possible bulk perturbations - in particular, we include 
SU(2) breaking terms. Such bulk perturbations are im- 
portant to analyze because they may qualitatively change 
the starting point for the boundary RG from which the 
A (g> N fixed point was derived. At the free fermion fixed 
point, the leading physical perturbations are marginal: 
they involve various combinations of bilinears in d<fi a 
and e ±l< ^° , where <f>J , a = 1, 2 are the two left/right 
moving bosonic fields in an appropriately normalized ba- 
sis. Some of these operators correspond to changes in 
Luttinger parameters as well as perturbations away from 
the equal velocities case. In the free fermion limit this 
does not open any bulk gaps, and we conjecture that this 
is the case even with interactions. Hence they likely do 
not destabilize the boundary A ® N fixed point. There 
is also the operator cos(v2#<r), which corresponds phys- 
ically to an attractive U Hubbard interaction inducing 
power-law superconducting long range order and the for- 
mation of a spin gap. At this bulk RG fixed point electron 
tunneling onto the Majorana mode becomes highly irrel- 
evant: the electrons are bound up in Cooper pairs. Of 
course, the entire lead has then already become super- 
conducting, so there is no meaningful way to discuss a 
zero bias anomaly due to Majorana modes in this case. 
The remaining operators break the SU(2) spin symmetry 
to a Z 2 Ising symmetry which in turn is spontaneously 
broken, resulting in a ferromagnetic bulk. This situation 
is dual to the one in which cos(v / 2#o-) condenses, so the 
A®N fixed point does not survive here either. However, 
for all physical SU(2) breaking perturbations, such as 
spin orbit coupling and magnetic field, the A ® N fixed 
point is stable. 



superconductor one can also make the normal lead either 
helical (effectively spinless) or spinful. We found that the 
transport signatures of such junctions are radically dif- 
ferent and allow a robust distinction between topological 
and trivial superconductors. 

Specifically, we found distinct low energy fixed points 
corresponding to perfect normal and perfect Andreev 
reflection, and determined their regimes of stability in 
the different junction archetypes. Crucially, for the 
most physically relevant case of weak repulsive interac- 
tions we found that the stability of the normal reflect- 
ing fixed point is equivalent to the absence of a Majo- 
rana zero mode. Conversely, the presence of a Majo- 
rana zero mode in this regime can be uniquely detected 
through a quantized conductance G — 2e 2 /h character- 
istic of the Andreev fixed point. We showed that such 
a quantized conductance is indeed a 'smoking-gun' sig- 
nature of topological superconductivity. It should be 
again stressed that this distinction is not necessarily cap- 
tured in the free fermion case — scattering theory in the 
non-interacting limit allows for similar conductances in a 
spinful/topological and spinful/ non-topological junction. 
Fortunately, repulsive interactions restore the sharp con- 
ductance dichotomy between these setups. 

Let us summarize the physical consequences of the 
renormalization group flows for each of the junction 
archetypes. 



Helical wire — topological superconductor 
junction 



Let us first discuss the helical/topological junction, 
which in the regime g > 1/2 should manifest a zero- 
bias anomaly associated with tunneling onto a Ma- 
jorana fermion. Typically captured in terms of free 

"""In 



fermion pH3oMi33i35] ( but see R e fe, 1551571) , we f oun d the 
zero-bias anomaly to be robust to rather strong repul- 
sive interactions and arbitrarily strong attractive inter- 
actions. To calculate the value of the zero-bias conduc- 
tance, we compute the imaginary time current-current 
correlation function 



VI. DISCUSSION 



H(r) = (I(r)I(O)), 



(89) 



In light of the intense current effort to realize Majorana 
fermions in condensed matter systems, it is important to 
understand the experimental signatures of these topolog- 
ically protected zero modes. In this work, we focused on 
tunnel junctions from normal Luttinger liquid leads to 
superconductors in order to probe the existence of Majo- 
rana zero modes in the latter. Our approach was based 
on universality, and the low energy fixed points we found 
govern a wide variety of experimental setups — not just 
junctions in the limit of weak tunneling. In particular, 
from a low energy universality point of view, it makes 
sense to consider four general junction archetypes, where 
in addition to the binary choice of topological or trivial 



which determines the conductance through the Kubo for- 
mula 



G 



h 



iS) 



(90) 



Since @/tt gives the total charge on the Luttinger liquid, 
the current is I — eO/n. Using the fixed-point bound- 
ary action of Eq. (36 1 (with v = 0) it is straightforward 
to show that H(iui) = e 2 g\uj\/ir, which apparently yields 
G = g (2e 2 /h). However, this result fails to take into ac- 
count the finite extent of the helical wire and the fact that 
it must ultimately be contacted by Fermi liquid leads [see 
Eq. ( 36 1] . This is a familiar problem in Luttinger liquid 
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FIG. 5: Sample dl/dV curve as a function of voltage V (in 
arbitrary units) for g = 0.95 in a helical nanowire/topological 
superconductor junction. The blue and red portions are the 
perturbative results around the normal and Andreev fixed 



points respectively, calculated according to Eq. (91 1 



theory, and can be resolved in various w ays by correctly 
modeling the wire together with the leads 58 59 . A partic- 
ularly convenient scheme is to continuously interpolate 
the Luttinger parameter between its interacting value of 
g and its free fermion value of 1. [One can implement 
this interpolation scheme with a spatially varying g(x) 
that takes on a value of g close to the junction and 1 
far away. In our boundary theory for the junction, this 
produces a frequency-dependent g(oj) that goes to one 
at zero frequency and g at high frequency.] This simple 
model already exhibits a crossover of the frequency de- 
pendent conductance G(ui) from g (2e 2 //i) for u > Vf/L 
to 2e 2 /h for cj < Vp/L (assuming this energy scale is 
smaller than the other relevant scales in the problem). 
Hence, even in interacting junctions, the zero-bias con- 
ductance is expected to be quantized at G = 2e /h. 

Universality also gives the line-shape of the finite bias 
conductance curve dl/dV as a function of V. The RG 
flow from the unstable normal reflecting fixed point to the 
stable Andreev one defines a crossover voltage V* , which 
is roughly the width of the zero bias anomaly. The entire 
dl / dV curve can in principle be determined from univer- 
sality. Although the general calculation is involved, it is 
easy to calculate dl/dV for V > V* and V < V* by 
performing perturbation theory around the normal and 
Andreev fixed points respectively. We obtain 



G 



(V/F*)-2(i-i/(2 fl )) ) v^V* 
2e 2 /h-(V/V*) 2 ^~ 1 \ V^V* 



(91) 



with the corresponding conductance line-shape illus- 
trated in Fig. p)J Equation (91 1 is valid down to voltage 
hvp/(eL), which we assume to be much smaller than V* . 
Extracting the same value of g from fitting experimental 
results to Eq. ((Ml for both high and low V would provide 
a non-trivial check of our results. 

In fact, there is some subtlety involved in the derivation 
of (91). Indeed, for V ^> V* , general scaling arguments 



G(V,t) = G 



t/t 



(V/Vo) 1 - 1 /(2s) 



(92) 



where G is a scaling function and to is defined at the 
cutoff scale Vq- G(x) is an even function of x, so gener- 
ically we exp ect G(x) ~ x 1 . However, for free fermions, 
(13 1 and (10 1 show that G(x) ~ x 4 , i.e. the coefficient 



of the quadratic term in G vanishes. We believe this to 
be highly non-generic, and expect that as soon as inter- 
actions are turned on, the coefficient becomes non-zero. 
We can also consider finite temperature, in which case 
we have a two-parameter scaling function 



G(F,T,t) = G(t/U 1/2 ,t/r 1/2 ) 



(93) 



Again, while in the free fermion case G{x) has a non- 
generic vanishing of the quadratic term in x, the lin- 
ear response conductance at high T has a non-vanishing 
quadratic term. That is, G with the first argument set 
at infinity has non-zero second derivative with respect 
to the second term. This qualitative distinction between 
zero temperature free Fermi and generic results again un- 
derscores the power of universality and the utility of our 
approach. 

Apart from the conductance, the fixed points studied 
here can be further distinguished by the behavior of the 
local density of states at the junction. Following Ref. 55 
the junction's local density states at frequency u> evalu- 
ated at each fixed point is given by 



PldosM 




(perfect normal reflection) 
(perfect Andreev reflection). 

(94) 
For the helical/topological junction, normal reflection is 
stable in the regime g < 1/2, and the local density of 
states vanishes as a power law in u). More interestingly, 
for 1/2 < g < 1, which is likely the most physically 
accessible regime for solid state systems, the Andreev 
fixed point is stable and results in a divergent density 
of states at low energies. This remarkable feature has 
also been predicted for a spinful Luttinger liquid-ordinary 
superconductor junction with perfect Andreev reflection 
boundary conditions (but see the discussion below), and 
may be observable in tunneling measurements 5 - 5 -! At- 
tractive interactions corresponding to g > 1 remove this 
divergence, despite the perfect Andreev reflection fixed 
point remaining stable there as well. 



B. Helical wire- 



-non-topological superconductor 
junction 



show only that 



In the case of a helical/non-topological junction, the 
normal reflecting fixed point is stable for g < 2. Thus, 
in the physical regime of weak repulsive interactions, the 



19 



h dl 0.12 
2^dV „,„ 



0.06 
0.04 



V 



V 



FIG. 6: Sample dl/dV curve as a function of voltage V 
(in arbitrary units) for g — 0.7 in a helical nanowire/non- 
topological superconductor junction. The red portion is cal- 
culated according to Eq. (951, while at sufficiently low volt- 



ages (blue curve) the Fermi liquid leads dominate and we have 
a crossover to g — 1 scaling. 



zero-bias tunneling conductance is predicted to vanish. 
We can again use perturbation theory to extract the finite 
bias conductance for small V. We obtain the power law 
vanishing form 



G~ (V/V*) 



2(2/ff-l) 



(95) 



SU(2) invariance) for simplicity, we obtain 



G 



(V/V*) s*-, V»V* 

2e 2 /h-(V/V*)^+ r , V<t:V* 



(97) 



The local density of states at the junction in the 
physically relevant regime of g p < 1 can be computed 
by evaluating the fermion two point function. For the 
spinful/non-topological junction, we must evaluate the 
two point function at the normal reflecting fixed point, 
resulting in 



PLDOS H ^ 1/(29p) - 1/2 , 



(98) 



which goes to zero at low energy. This is a generic fea- 
ture of the g p < 1 regime of the spinful/non-topological 
junction, resulting from the instability of the perfect An- 
dreev reflecting fixed point. On the other hand, for the 
spinful/topological junction, the two point function must 
be evaluated at the Andreev fixed point, at least for 
g p > 1/3, resulting in 



PLDOSH-W 9 ' 72 - 172 - 



(99) 



The local density of states thus diverges for weakly re- 
pulsive interactions in this case. 



We expect this form to be valid down to voltages eV ~ 
hvp/L, below which point the finite length of the inter- 
acting Luttinger liquid becomes important. Indeed, if 
we model the helical wire as a Luttinger liquid whose 
Luttinger parameter changes continuously from its inter- 
acting value of g to the free Fermi value 1 at x ~ L, then 
for eV ~ hvp/L we are mostly sensitive to the long dis- 
tance free Fermi part of the lead, and thus we expect the 



conductance to cross over to the behavior given by ( 95 ) 
with g = 1 (see Fig.|6|. 

The stability of the normal reflecting fixed point in the 
physical regime 1/2 < g < 1 also yields a different local 
density of states for the helical/non-topological junction. 
Indeed, Eq. (pi| shows that 



/OldosM ~ wl/s * 



(96) 



vanishes at zero energy. 



C. Spinful junctions 

The analysis of the spinful junctions proceeds analo- 
gously to that of the helical junctions. In the case of a 
topological superconductor we once again have a zero- 
bias anomaly, and we can compute the shape of the con- 
ductance curve for high and low V from the scaling di- 
mensions of the leading perturbations to the normal and 
Andreev fixed points respectively. Using the results of 
the previous section, and assuming g a = 1 (i.e., spin 



D. Future directions 

Although we have shown that perfectly quantized con- 
ductance is a universal property of the low-energy/long- 
distance limit in topological junctions, a relevant issue for 
experiments is how closely one can approach this limit in 
practice. In particular, in our analysis we have assumed a 
semi-infinite Luttinger liquid lead, and an arbitrarily long 
topological superconductor. In a physical setup, neither 
assumption is valid. In particular, a finite topological su- 
perconductor of length Lgc wm have additional couplings 
6 and t' (c.f. the discussion at the end of Sec. Ill A) of 
magnitude ~ exp(— £sc/£)j where £ is the induced co- 
herence length. These couplings are relevant (for g < 2 
in the helical case) and ultimately drive the system to 
the normal reflecting fixed point, with zero conductance. 
However, if they are small to begin with, the crossover 
will occur at low energies and an intermediate regime 
with enhanced conductance is to be expected. 

There exist many directions for further investigation. 
One has to do with the nature of the crossover of the 
quantized conductance from G = 2e 2 /h to G = g(2e 2 /h) 
as a function of the driving frequency uj, which is ex- 
pected to occur for u> ~ wf/^ll, where Lll is the length 
of the Luttinger liquid lead. Extracting the full depen- 
dence on both to and the temperature T is subtle and 
may require use of the Keldysh formalism. Another is to 
study the fate of the Majorana zero mode for g < 1/2 in 
the helical/topological junction. Here we expect the Ma- 
jorana mode to be hybridized with the degrees of freedom 
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in the wire but nevertheless to remain power-law local- 
ized near the junction. It would be interesting to study 
the form of this Majorana zero-mode using DMRG sim- 
ulations. Yet another direction would be to extend the 
present analysis to the case of the multi-channel case, 
which may be an experimentally relevant regime. An- 
other potentially experimentally relevant issue is the role 
of disorder in the Luttingcr liquid lead, and its effect on 
the conductance. Finally, it would be interesting to in- 
vestigate multi-terminal junctions and possible universal 
signatures here. 
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Appendix A: Derivation of the effective action for $ 

We start with the action 



S = drl H(<f>, 



dx d T 8d x 



with 



r°° ,, 

H = J dx—^id^Y+g-^dJ) 2 



(Al) 



(A2) 



The normal reflection fixed point pins 9(x = 0). Consider 
the right and left moving fields: 



4> R = dp + 9/g, 4>l = 4>- 0/g 

The boundary condition 6(x = 0) = gives 

4> L (x = 0) = <f> R (x = 0). 

We define a new field <p(x) as 

<j>{x) — <j>n{x), x > 
4>{-x) = (J>l(-x), x < 

Note that 0(0) = <j> R (0) = <j> L (0) = 0(0). We have 



So = g / dxdr 



iiM? + i& 



(A3) 

(A4) 



(A5) 
(A6) 
(A7) 



(A8) 



We introduce a variable ( I > (t) and a Lagrange multi- 
plier A as 

S A = i f drA(r)($(r) - 4>(x = 0, t)) (A9) 

We now integrate over the modes 0fc )W . Note that 
dk diu 



Sn = 



(27T) 



7 G (k,u)(/)k,w<f>-k,-w 



with 



G -1 (ife,w) 



gv 

-in 



2v 4w 2 



(A10) 



(All) 



Performing the gaussian integral over 4>k,u yields 

dk diu 1 



duj 
2tt j 



(2tt) 2 4 
Performing the integral over k gives 



G(fc,w)A ClJ A_ u; (A12) 



dto 

2^ 



iK$-u> 



2g\u 



K^-w 



(A13) 



Finally, performing the functional integral over A w yields 



S\ 



2tt 2tt 



(A14) 



as desired. 



Appendix B: Derivation of the path integral at the 
perfect normal reflection fixed point 

Here we will derive the path integral representation for 
the partition function in Eq. ( 30 1 , beginning from the full 



action that includes bulk degrees of freedom in the helical 
Luttinger liquid. This is done by inserting resolutions of 
the identity between intermediate imaginary time steps, 
using the dual orthonormal bases 

1 = 5>^>< S ^|=5>*0)( S *0|. (Bl) 

Here s x,z denote the eigenvalues (= ±1) of the corre- 
sponding Pauli operators a x,z , and <f>,6 are shorthand 
for eigenvalues of <fi(x),6(x),0 < x < L. To obtain the 
partition function we need the overlaps between these 
states. Using 

|«« = ±1> = -^(|«* = 1)± !«* = -!)) (B2) 



we obtain 



(s x \s Z ) = -^exp 



(1 -s x ) (l-s z ) 
in 



(B3) 
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which together with the fact that <j>(x) and d x 0{x)/'K are 
dual variables yields (up to unimportant constants) 



Appendix C: Solution of a non-interacting helical 
wire coupled to a single Majorana mode 



>s x \9s 
= exp 



(l-s x ) (1 - s z ) i 



dx 4>d x 9 



(B4) 



2 2 

It follows that the partition function can be expressed as 



Z 



v^ve ]T J2 



-Se+iSb 



(B5) 



{s=(T)e±l}{s^(T)6±l} 

with Se the usual Euclidean action 
"1 



Se = — dxdr 

2tt 



2t / drcos 
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-{d x 6) 2 + g(d x ^y 



(X = 0) + 7T 



1 



(B6) 



and Sb the Berry phase contribution 
5s = 7 f dr{l - s z )d T s x + - f dxdrd T (j)d x e. (B7) 



Suppose that we sum over s x in Eq. ( B5 ) . To do so it 



is convenient to first rewrite the path integral exchanging 
4>(x, t) for the new variable 



(x,t) 



(x,r) + -(l-s x (r)) 



(B8) 



Integration by parts shows that s x then appears only in 
the following contribution to the action: 



dr^ld r 



<x = 0) + - (1 - s z ) 



(B9) 



Summing over s x , at each discrete time step, gives a 
vanishing contribution to the partition function unless 
d T [~9(x = 0) +tt/2(1- s z )} = modulo 2vr at each 
t. Hence the exponential of 2iri times this quantity 
is conserved. Converting to operator notation, this is 
<t z (t) exp [i6(0, t)] — which is just the fermionic parity of 
the superconductor plus Luttinger liquid system. 



Let us now instead sum over s z in Eq. (B5 I. Using Eq 



(B7l, it follows that s x is a conserved quantity, so the 



path integral sums over only two imaginary time config- 
urations of s x : s x = +1 for all r and s x = — 1 for all 
t. Integrating out everything but $ = <j)(x = 0) yields 
the action S noTma \ + 5St derived above [see Eqs. (20l 28 1] . 
The partition function then becomes 



J „X L 



-Snormal p ~2t f (ITS* COS $ 



(BIO) 



s x =±l 



which indeed recovers Eq. ( 30 I 



In this Appendix we will sketch the solution of the 
following Hamiltonian, 



H = 



dx 



iv F ^d x ip + ^=7i(V> f - ip)S(x) 



(CI) 



which describes a non-interacting helical wire of length L 
coupled to a single Majorana mode 71 at one end. This 
Majorana's 'partner' 72 is assumed to decouple entirely 
from both the wire and 71. At the other end of the 
wire we will impose perfect normal reflecting boundary 
conditions, requiring 



ip(x — L) = ip(x = —L). 



(C2) 



Our goal will be to find the low-energy wavefunctions of 
H in the limit E <C t and expand ip and 71 in terms of 
the corresponding modes. 

As a first step we write H in terms of W(x) — 
[^(x)i>(x)Pf], where / = (71 + 172) /2: 



H = 



H = 



/ dx&Wb. 






—ivpdx 

—ivFd x 
-frix) ^5{x) 
-^5{x) j-6{x) 











(C3) 



(C4) 



The wavefunctions <&e(x) with energy E can be immedi- 
ately written for x ^ as 



$ E (x > 0) = 


r *f* > " 

e »^o| 




d E 



,$e(x<0) 



e "^ 

ce 
d E 



(C5) 



The elements above are constrained by normalization, 



the boundary condition of Eq. ( C2 ) which requires 



e "-f a^ 






"*■ a f 



°E' 



(C6) 



and the following relations needed to satisfy the Hamil- 
tonian at x = 0. 





Ec E 



= (a> - a<) + (&> - b<) 

t 



iv F (a E 



a E ) + -pji c E + d E ) 



(C7) 



7t v 



= £d B = — (6>-a<). 



Equations (|C6| and (|C7|) admit a non-trivial solution 

. (2EL\ 



provided the energies satisfy 
Ev F 



= 



f 2 



[2EL\ 
cos — 1 



{-^ 



(08) 
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For E <t the energies are well-approximated by 

nirvF 



E n 



2L 



(C9) 



for integer n; the associated wavefunction components 



can be found to leading order in E n /t from Eqs. (C6) 
and dCTl. 



One can deduce by inspection that one of the zero- 
energy wavefunctions supported by H is 



ft, 





" " 


1 





/2 


i 




— i 



(CIO) 



which simply corresponds to the decoupled Majorana 
mode 72. The wavefunctions carrying energy E n with 
n even are given by 



$/ 



(n even), 



(Cll) 



where JV n is the normalization. Note that §e q corre- 
sponds to the zero-energy Majorana mode which is ab- 
sorbed into the helical wire due to the tunneling t. Fi- 
nally, the n odd wavefunctions read 



$B.. = 



isgn(x)e v f 

- E n 

— isgn(a;)e "* 

t 



(n odd). (C12) 



With these wavefunctions in hand one can now expand 
ip(x) and 71 in terms of low-energy modes for the system: 



ijj(x) 



n even 



7i = P+f 




(C13) 



Here r+ = r_„ and TjJ = T'_ n respectively create energy 
E n excitations with n even and odd. Equations (C13) 



encode two important relations. First, it follows that at 
low energies 

^(x = 0+) = iP R (x = 0) = V f (x = 0~) = ^ L {x = 0), 

(C14) 
which is the familiar perfect Andreev reflection bound- 
ary condition induced by the coupling to 71. Second, 
71 and i[>, ift are not independent at low energies; using 
Eq. ( C14| their relation can be expressed in the following 
symmetric form 

71 - i ^|[v t (o + )-^(o + )-V- t (o-) + V'(o-)] 

= i-%[^(0) ~ MO) - 4(0) + MO)]- (C15) 



Appendix D: Solution of an interacting helical wire 
coupled to a single Majorana mode 

We now treat the interacting helical wire. We first 
redo the calculation for the non-interacting case directly 
in terms of the bosonic modes, and then generalize to 
5 ^ 1. It is useful to bosonize both the Luttinger liquid 
and the two-level system formed by 71 , 72 , as done in 
Section HTU 



7i 

72 



a" 



(Dl) 
(D2) 



Also, because the Jordan- Wigner string goes to the left 
in ( |26[ ), the bosonized form of the continuum Fermi 
fields includes an extra factor of a z , which is simply the 
Fermionic parity of the 71, 72 system. The tunneling term 
is H t = 2ta x cos((f>(x = 0)), but for convenience in this 
Appendix we shift the phase of cf>(x) by ir/2, so that the 
tunneling term becomes 



H t = 2ta x sin((t)(x = 0)). 



(D3) 



From now on we work in the low energy Hilbert space 
"Hiow C H where <fi(x = 0) is pinned at ±7r/2, i.e. we 
assume we are exactly at the Andreev fixed point. The 
state of the spin-1/2 representing the topological super- 
conductor is then completely determined by 4>(x = 0), so 
that we need only retain the Luttinger liquid degrees of 
freedom to describe all the states and operators in T-L\ ow . 

In particular, the operator 72 = a x , representing the 
decoupled Majorana mode, is given by sm(<j)(x = 0)). 
Now, if the fermionic Hamiltonian were quadratic, we 
would necessarily have a partner Majorana mode for 72, 
i.e. another fermionic operator which commuted with the 
Hamiltonian and squared to 1. Indeed, this operator is a 
zero momentum "plane- wave" solution of the Bogoliubov- 
de Gennes equation (see Eq. Cll). What about the 
interacting case? In this section we derive a bosonized 
expression for the partner Majorana mode in such an 
interacting helical wire, for all g > 1/2. Furthermore, 
we check that this expression is correct perturbatively to 
leading order in g — 1. 

We begin by redoing the calculation for g = 1 in 
the bosonic framework, and then generalize to g =/= 1. 
The normal reflecting boundary condition at 6(L) = is 



equivalent to 



94>(L) 



= 0. A general <fi field configuration 



can thus be expanded as 



(x) = ±tt/2 + 2_. 4>n sin 



(2n + l)irx 



The Hamiltonian for g = 1 reads 



H = 

Expanding 
d 



2tt 



dx 



{d x <j>? 



2L 



d<j)(x) 



L , . f(2n + l)Trx\ d 

dx sin — - ; . . 

\ 2L J dcb(x) 



(D4) 



(D5) 



(D6) 



T>> 



and inverting, we obtain 

n = ^ E 



ttvf 



n>0 



a \ 2 2 



(D7) 



with w n = (2n + l)/4. We define creation and annihila- 
tion operators 



1 







2u>, 

LO, 



(o„ + 4) 



and expand out the field operators 4>(x) and 
6'{x) = -7 



dy 



84>{y) 



(D8) 
(D9) 

(D10) 



Here 9'(x) = 9(x) — 9(L) sends the Jordan- Wigner string 
to the right, and thus preserves the boundary condition 
0(0) = ±7r/2. We obtain: 



e\ x ) 



2 2^ 






/2u} r , 



E 






/2uj,-, 



H 



(Dll) 

(D12) 



where k n = 2iru> n /L. We can exchange 8 for 9' at 
the expense of introducing an extra factor of Pll = 
exp(#(0) — 0(L)), so that the bosonization now becomes 



IpR/L 



1 - 26l 1 6_ 1 ) e^ ±e > 
P LL e^ ±e 'K 



(D13) 



For the sake of efficiency, we now change notation 
(tl!ii(x),ipL(x)) -^ ( r ( a; )j^( :c ))- For a quadratic Fermionic 
Hamiltonian (g=l), one can easily verify that 



<h 



i 
~4L 



dx (r - r 1 + I - £ f ) 



(D14) 



is the partner "plane- wave" Majorana mode for 72. 8\ is 
the result of 71 leaking into the Luttinger liquid via the 
coupling t. with (D14l valid in the low energy Andreev 



limit. (D14l bosonizes to 



61 = 



Pll 

2L 



dx [sin(</> + $') + sin(0 - 6')] . (D15) 



u 



It is instructive to expand out (D15l in oscillator modes 



and verify explicitly that the bosonized expression com- 
mutes with the bosonized Hamiltonian, anti-commutcs 
with 82, and squares to 1. We represent 8\ as 



*i = 



<5f 



(D16) 



with respect to the decomposition of T~L\ ow into the two 
(f)(0) = ±? sectors. Note that S\ is purely off-diagonal 



because the factor of Pll in (D15l exchanges the two 



sectors. This immediately shows that 5\ anti-commutes 
with 



5o = 



I 
-/ 



(D17) 



The bosonized expressions for both 8^ and 6 1 + are sim- 



ply (D15l with the factor of Pll stripped off, but it is 



useful to re-write them as 



dx 



cos((</> - §) + 9') + cos((0 - §) - 9') 



6i + = - dx 
Jo 



2L 

<- L cos((0+f) + e')+cos((0+f)-0') 



2L 



This form removes the constant ±^ from <j) and in par- 
ticular shows that only terms with an even number of 
creation/annihilation operators appear with non-zero co- 
efficients. Furthermore, the coefficient of a potentially 
energy- violating term would be the real part of an oscil- 
latory integral which is purely imaginary; hence 8\ com- 
mutes with the Hamiltonian. Similarly one can compute 
8\ directly and see that it equals 1, though it is instruc- 
tive to check this explicitly on some low energy subspaces. 
The lowest lying states are 



V2 



Krio) 



{-^(4) 3 |0),ai|0)} 



E = 





E = 


TTVp 


E = 


TTVp 

~2L 


E = 


Sttvf 



4L 



(D18) 



One can explicitly compute that on the lowest 3 sub- 
spaces 5± acts as +1,-1,-1 respectively, whereas on 

4L 



(D19) 



the E = 3 Tr F subspace it acts as 




a non-trivial matrix that squares to +1. 
When g 7^ 1, we have 



H = 



■kvf 



E 

n>0 




(gu n ) 



.(D20) 



We can still expand </> in modes as in ( D4 ) , but this time 

1 



V^g^n 



(o„ + <) 







gu>n 



(D21) 
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To get rid of unwanted g's we define 

<^0z)=±|+Vs/ dydyfty). (D22) 

4> 9 (x) is thus diagonal with respect to the <^>(0) = ±7r/2 
sector decomposition, with the choice of sign in (D22) 



corresponding to the choice of sector. Essentially, we 
have used the fact that (f)(0) is pinned to give well-defined 
meaning to the operator y/g<f)(x). Similarly we can use 
the fact that 9'(L) — to define 



o/y 



(x) 



V9 Jx 



dyd x 6'(y) 



(D23) 



The virtue of the operators <fr 9 (x) , 6' ' 9 (x) is that their 
mode expansions are i denti cal to thos e of cf>(x),6'(x) in 
the free case, given by (Dll I and (D12). Thus Sf, defined 

by 



5 1 = -^r I dx [ smi 



+ i 



sM4> 9 - e' 9 ) 



2£ jo 

(D24) 

has all the necessary properties: it anti-commutes with 

82, commutes with H, and squares to 1. It can also be 

written in terms of the fermions. Let us see how that 

works explicitly for small e = g—1. We expand to leading 

order: 



9 {x)^(t>{x) + 



'(as) « 9'(x) + 



dyd v 4>{y) 



dyd y e'{y), 



(D25) 
(D26) 



so that 

shi{(/) 9 + 6' 9 )^sin((t) + 9') 

dyd y (/)(y) 



dyd y 6'{y) . 
(D27) 



We now rewrite (D27) in terms of fermions using 

d x <j> = 7T (rV - lU) (D28) 

dj = it (rV + lU) , (D29) 

to obtain: 



P LL sin(^(x) + 6' 9 (x)) = ~-(r(x) - r\x)) 



7re , 



--^(r(x)+r^(x)) 



dy(r i r + tf£)-2 / dylH 
J a 



(D30) 



A subtle point is that (D28l has to be normal ordered 



with respect to the ground state of the free fermions with 
the tunneling term included; in practice this means re- 
placing r'r — > (r^r + rr^)/2, etc. With a similar ex- 
pression for Pll sin(</> 9 (x) — 6' 9 (x)) we finally obtain the 
interacting Majorana mode 



S 9 = S x + 6[ 



(D31) 



with 



S 'i = ^T I dx(r + r^ -£-£ r ) f dy{r^r + tfl) 
8L Jo Jo 

dx [ dy[{t{x)+fi{x))r^r{y) 

(r(x)+r\x))tU(y)]. (D32) 



7re 

4L 



L 
1 
^0 



Likewise we can also expand the Hamiltonian in 
Fermions: 



H = H Q + H' 



(D33) 



with 



Ho 



and 



IVp 



dx (-r^d x r - rd x r ] + tfd x £ + £d x fi) 



H' = -2i:ev F / dxr^rtH. 
Jo 



(D34) 



(D35) 



Note that the Hq we have differs from the usual one by 
a boundary term; again this boundary term effectively 
takes into account the tunneling at x = 0. We now need 
to check [H, Sf] = 0, which at leading order reduces to 
checking that 



[H ,5' 1 ] = [6 1 ,H'] 



(D36) 



( D36 ) can be explicitly verified with some algebra, show- 



ing that ( D32 1 is the correct leading order interacting 



correction to the Majorana mode S\ 
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